Draft version February 2, 2008 

Preprint typeset using I^T^jX style cmulateapj v. 10/09/06 



HST MORPHOLOGIES OF LOCAL LYMAN BREAK GALAXY ANALOGS I: 
EVIDENCE FOR STARBURSTS TRIGGERED BY MERGING 

Roderik A. Overzier 1,2 , Timothy M. Heckman 1 , Guinevere Kauffmann 2 , Mark Seibert', R. Michael Rich 4 , 
Antara Basu-Zych 5 , Jennifer Lotz 6 , Alessandra Aloisi 7 , Stephane Charlot 8 , C. Hoopes 1 , D. Christopher 

Martin 9 , David Schiminovich 5 , Barry Madore 3 

Draft version February 2, 2008 

ABSTRACT 

Heckman et al. (2005) used the Galaxy Evolution Explorer (GALEX) UV imaging survey to show 
that there exists a rare population of nearby compact UV-luminous galaxies (UVLGs) that closely 
resembles high redshift Lyman break galaxies (LBGs). We present HST images in the UV, optical, 
and Ha, and resimulate them at the depth and resolution of the GOODS/UDF fields to show that the 
morphologies of UVLGs are also similar to those of LBGs. Our sample of 8 LBG analogs thus provides 
detailed insight into the connection between star formation and LBG morphology. Faint tidal features 
or companions can be seen in all of the rest-frame optical images, suggesting that the starbursts are 
the result of a merger or interaction. The UV/optical light is dominated by unresolved (^100-300 
pc) super starburst regions (SSBs). A detailed comparison with the galaxies Haro 11 and VV 114 at 
z = 0.02 indicates that the SSBs themselves consist of diffuse stars and (super) star clusters. The 
structural features revealed by the new HST images occur on very small physical scales and are thus 
not detectable in images of high redshift LBGs, except in a few cases where they are magnified by 
gravitational lensing. We propose, therefore, that LBGs are mergers of gas-rich, relatively low-mass 
(M* ~ 10 10 Mq) systems, and that the mergers trigger the formation of SSBs. If galaxies at high 
redshifts are dominated by SSBs, then the faint end slope of the luminosity function is predicted to 
have slope a ~ 2. Our results are the most direct confirmation to date of models that predict that 
the main mode of star formation in the early universe was highly collisional. 

Subject headings: cosmology: observations - early universe - galaxies: high-redshift - galaxies: star- 
burst 



1. INTRODUCTION 

How did galaxies form? Ultimately, this simple ques- 
tion captures most, if not all of the most widely pursued 
research topics in cosmology. Surveys of galaxies at dif- 
ferent redshifts are being used to study which galaxies 
contain most of the stellar mass at a given epoch, and 
which galaxies are undergoing the strongest evolution. 
These measurements serve to constrain models of struc- 
ture formation. The main formative processes associated 
with typical galaxies (e.g. star formation, merging, and 
feedback) were largely completed within the first half of 
the Hubble time, making their study a very challenging 
one. One of the best probes for studying star forma- 
tion in the early universe is provided by galaxies that are 
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luminous in the rest-frame UV due to intense star for- 
mation. The most luminous of these, the Lyman break 
galaxies (LBGs), are easily detected at z = 2 — 6 in deep 
pencil beam surveys from the ground and with the Hub- 
ble Space Telescope CHSTl (e.g.lSteidel fe Hamiltonlll993l: 
I Steidel et al.lll999HShaplev et al.l|2001l: IGiavalisco et al.l 
20041: lAdelbereer et alJ 12004 Bouwens et alJ 120061 : 
Yoshida et a 1.1 12006D . The cosmic star formation rate 
(SFR) reached a maximum at z ~ 1 — 3 and de- 
creased dramatically toward z — (e.g. iLillv et aTl 
119961 : iSchiminovich et afll2005h . There is only a mod- 
est decrease in the UV luminos i ty de n sity from z ~ 
3 to z ~ 6 (e.g. lOuchi et al.l 120041 : iBouwens et~aT1 
2006), indicating that LBGs represent a major phase 
in the early stages of galaxy formation and evolu- 
tion. Based on their clustering and number statistics, 
LBGs are the precursors of present day massive galax- 
ies undergoing a phase of intense star formation, and 
a large fraction of LBGs may merge to form ellipti- 
cal galax ies that are situated in present-day groups and 
cluste rs dGovernato et all 120011: IGiavalisco fc Dickinson! 
1 20011: iMoustakas fc Sqmervillel l2002t lOuchi et all 12004 
lAdelbereer et al.ll2005f) .~ 

Studies of the structure and sizes of LBGs indicate 
that high reds hift galaxies are compact (^0.1 // -0.3 // , 
or ^ 1 -2.5 kpc: | Giavalisco et al.lll996al lLowenthal et al.l 
119971 iFerguson et all l2004t IBouwens et all l2004h . and 
that large (>0.4") low surface brightness galaxies are 
rare. The morphologies are best characterized as being 
dominated by one to several UV-bright knots embedded 
in diffuse emission, with little difference between the rest- 
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frame UV a nd optical light distributions (Pa povich et al.l 
20011 1200^ . The morphologies of LBGs are unlike those 
of nearby galaxies on the Hubble sequence. This is true 
even if the images of the nearby objects are degraded 
so that they have similar resolution and depth to those 
of LBGs. LBG morphologies are more similar to those 
of ultraluminous infrared galaxies (ULIRGs), which oc- 
cupy a region of the concentration/asymmetry plane that 
lies betwee n double-nucleated mergers, spheroids and 
disks (e.g I Abraham et all Il996t [Conselice et al.l 12004 
lLotz et all 12004 120061 : lLaw et al.l l2007al and references 
therein). However, the interpretation of the UV mor- 
phologie s of high redshift populations remains problem- 
atic (see lLaw et al.ll2007al) . because it is unclear whether 
the irregular morphologies are a consequence of the merg- 
ing of gas-rich galaxies, or whether one is just seeing 
patchily distributed star formation within a single sys- 
tem. The answer to this question and how it evolves 
with redshift is essential for understanding how galaxies 
formed. 

A detailed study of LBGs is limited by two factors. 
First, their distances render them small and faint, and, 
second, star forming galaxies at high redshift are system- 
atically different from starburst galaxies in the nearby 
universe. Local starbursts (SBs) show strong correla- 
tions among their basic properties that are different from 
LBGs. In particular, there is a systematic increase in the 
amount of dust obscuration {LfirI Lfuv) with increas- 
ing SFR, corresponding to a systematic increase in metal- 
licity with mass. LBGs are substantially less obscured 
than local SBs of similar SFR or mass (e.g. iReddv et al.l 
2006). The relatively low extinction of LBGs is also a 
characteristic of blue compact dwarfs (BCDs, alsoor HII 
galaxies) , but the latter have SFRs that are typically two 
orders of magnitu de smaller than LBGs and they have 
lower masses (e.g. [F anclh ct al. 1988; Tellcs & Tc rlevichl 
[19971 iHunter fc Hoffman! Il999l : IGil de Paz et al.l l2003h . 
Altho ugh BCDs may re s emble the fa i ntest LBGs 
(e.g. iMeurer et aT] 119951 : iNoeske et all l2006t) . and 
LIRG S and ULIRGS may resemble the dustiest LBGs 
(e.g. lAdelbereer k. Steidel 120001: iGoldader et all 120021: 
Daddi et alj 120051: IChapman et alj l2005t iHuang et all 
20051: Ivan Dokkum et al.ll2006[ L typical local SBs are dif- 
ferent from typical LBGs. LBGs have systematically 
lower metallicities and higher gas-mass fractions than lo- 
cal SBs with the same SFRs. The typical mechanisms 
that trigger starbursts might also be very different at 
high redshift compared to low redshift. 

Clearly, if we could find relatively nearby starburst 
galaxies whose properties are the best match to LBGs, 
we could study for the first time at a much higher 
physical resolution how the vigorous star formation and 
morphologies of typical high redshift galaxies are re- 
lated. However, local LBG analogs are currently very 
rare; as i ndicated by the Ga laxy Evolution Explorer 
(GALEX; iMartin et all I2005D imaging survey in the 
near (NUV, A~2250A) and far (FUV, A~1550A) UV, 
the co-moving number density of galaxies with FUV 
luminosities similar to LBGs has declined by a fac- 
tor of several hundr ed b etween z = 3 and z = 0. 
iHeckman et al.l (|2005l ) and iHoopes et alJ (|2007[ ) (hence- 
forward called "H05" & "H07") identified a rare pop- 
ulation of low-redshift (z < 0.3) galaxies with proper- 



ties remarkably similar to those of LBGs by matching 
sources in the GALEX a ll-sky survey with the Sloan Dig- 
ital Sky Survey (SDSS; I York et al.) |2000D spectroscopic 
sample, and selecting according to 2 criteria specifically 
desig ned to match the typical UV properti es of LBGs 
(e.g. ISteidel et al.l 119961 : IShaplev et al.ll2001h : L FUV > 
I fuv > 10 9 L Q kpc 2 , where L FUV is the 
FUV luminosity, XP\, and I fuv is the mean FUV sur- 
face brightness interior to the SDSS it-band half-light ra- 
dius (Ifuv — \L fuv I ' 7rr e,«) ■ By further limiting the 
sample to z < 0.3 and excluding broad line active galac- 
tic nuclei (AGN), this resulted in a sample of 31 "super- 
compact UV-luminous galaxies" (UVLGs). Interestingly, 
once selected according to these criteria, the UVLG sam- 
ple proved to be similar to typical LBGs in all other mea- 
surable properties (see H05 and H07). Further d etails on 
the sample are given bv lBasu-Zvch et all (|2007f ). who in- 
vestigated the radio to UV continuum properties of the 
sample and found that they follow the radio-far infrared 
correlation of normal star-forming galaxies, but gener- 
ally have less dust than other star forming galaxies with 
such high specific SFRs. In this respect, they are simi- 
lar to LBGs. This suggests that the process that causes 
star formation in the supercompact UVLGs differs from 
other local star forming galaxies, but may be similar to 
LBGs. Our sample of UVLGs ("local LBG analogs") is 
therefore highly valuable for understanding the nature of 
Lyman break galaxies at high redshifts 10 . 

In this paper (Paper I) we present the first set of high 
resolution HST images of 8 UVLGs and discuss their 
morphologies. In a subsequent paper (Paper II) we will 
carry out a more detailed study using a large HST data 
set to be observed in Cycle 16. The structure of this 
paper is as follows. In Sect. 2 we describe the new 
observations with HST. In Sections 3 and 4 we present 
the data, investigate UVLG morphologies both qualita- 
tively and quantitatively, resimulate our data at higher 
redshift, and compare with morphologies of star-forming 
galaxies and LBGs as well as two of the most nearby 
LBG-like galaxies known (Haro 11 and VV 114). We 
discuss the implications of our results for understanding 
the nature of LBGs at high redshifts in Section 5, fol- 
lowed by a summary of the results (Section 6). We use a 
cosmology [fi M) fi Aj /i] = [0.27,0.73,0.73] with H = lOO/i 
km s^ 1 Mpc" 1 . 

2. DATA 

2.1. HST observations 

We have observed 8 of the nearest (0.091 < z < 0.204) 
and brightest local LBG analogs of H05 & H07 with 
HST in Cycle 15. To date, 7 UVLGs have been ob- 
served with the Advanced Camera for Surveys (ACS) 
High Resolution Camera (HRC) in the filter F330W, and 
with the Wide Field Channel (WFC) through the filter 
F850LP. In addition, one UVLG has been observed with 
the Solar Blind Channel (SBC) in the filter F150LP, and 
with the Wide Field and Planetary Camera 2 (WFPC2) 
through F606W, given the new constraints following the 
failure of ACS during the course of 2007. Ramp fil- 
ter images with a central wavelength equal to that of 

10 See the Appendix for a detailed analysis that demonstrates 
that the interpretation and conclusions presented in the recent pa- 
per by Scarpa et al. (2007) are in fact in error. 
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redshifted Ha were also obtained for all but the latter 
source. Each filter probes star formation and morphol- 
ogy at a different characteristic timescale. The F330W 
filter (central wavelength A C «3334A with an effective 
width AA ~ 548A) is the most sensitive to star for- 
matio n over the past ~100 Myr (jLeitherer fc Heckmanl 
Il995f ) and provide the closest match to the GALEX 
NUV filter (A C «2315A, AA - 7301). The Ff50LP fil- 
ter (A c f«1614A, AA ~234A), available for one of the 
UVLGs, provides an excellent match to the GALEX 
FUV filter (A C «1530A, AA - 2551). The F606W 
(A c «6001 A) and F850LP (A C «9170A) images are well- 
matched to the SDSS images in the r- and i-band and 
probe older (>Gyr) stellar populations. The ACS ramp 
filters (FR647M with a width of 207A or FR782N with 
a width of 52A, depending on the redshift) probe red- 
shifted Ha from HH regions tracing the youngest stellar 
population (O stars with lifetimes <10 Myr), and allow 
us to search for evidence of galaxy- wide outflows of ion- 
ized gas. 

The targets were observed for one orbit per filter with 
ACS, and two orbits per filter with WFPC2. The F330W 
image of SDSS J092600.41+442736.1 failed to execute 
due to a guide star problem. The Ff50LP, F330W, 
ramp filter, F606W, and F850LP images consisted of 
3, 3, 2, 6 and 3 exposures, respectively, to facilitate 
the removal of cos mic rays. The images we re combined 
using Multidrizzle ()Koekemoer et al.l l2002t ) . producing 
registered, cosmic-ray free, geometrically corrected im- 
ages. The SBC and HRC images have a plate scale of 
0"025 pixel- 1 and a resolution of ~tf'07 (FWHM). The 
WFPC2/PC images have a plate scale of 0'.'046 pixel" 1 
with a resolution of «0'.'ll (FWHM). The WFC images 
have a scale of 0'.'05 pixel -1 with a resolution of «0'/12 
(FWHM). Magnitudes i n the AB system were mea- 
sured using SExtractor (|Bertin fc Arnoutsl [T996) from 
mAB = — 2.51og 10 (counts/T e2: p) + ZPT, where the ze- 
ropoints (ZPT) are 22.448, 24.085, 23.004, and 24.862 
mag for F150LP, F330W, F606W and F850LP, resp. All 
magnitudes were corrected for G a lactic extinction using 
the dust maps of ISchlegel et al.l (|1998| ). The reader is 
referred to Table [1] for a log of the observations, Table [2] 
for SFRs and stellar masses of the objects in our sample, 
and Table [3] for the main photometric data from GALEX 
and ACS. 

2.2. Continuum subtracted Ha images 

The Ha ramp filter images were continuum-subtracted 
as follows. First we scaled the ACS zsso image to an 
artificial ^775 continuum image, by determining the ratio 
of the flux of the 1775 continuum at the wavelength of 
redshifted Ha (measured from the SDSS fiber spectrum 
on opposite sides of the Ha line) to the total flux in the 
Z850 image. The total flux was measured within an aper- 
ture similar in size to the 3" diameter SDSS fiber aper- 
ture. We then created a "continuum-free" Ha image by 
subtracting the artificial 1775 continuum image from the 
ACS ramp filter image. The resulting Ha image should 
be a very good approximation, under the (very reason- 
able) assumption that the morphology of the continuum 
does not change between 1775 and zsso- I n the case of 
object SDSS J080844.26+394852.4 the subtraction was 
problematic due to the fact that it relied on subtraction 



of two bright point sources in z§50 and Ha. 

2.3. Spitzer observations 

Infrared photometry was obtained with the Infrared Ar- 
ray Camera (IRAC) and the Multi-band Imaging Pho- 
tometer for Spitzer (MIPS) aboard the Spitzer Space 
Telescope (PI: C. Hoopes, #20390). The total integra- 
tion time with IRAC at 3.6, 4.5, 5.8 and 8.0 /xm was 
60 s for each channel. The total integration time with 
MIPS at 24 and 70 /xm was 42 and 60 s, resp. We used 
a minimum 5 step dither pattern for removal of cosmic 
rays. The post-basic calibrated data (BCD) delivered 
by the pipeline were used to measure the integrated flux 
densities. The infrared photometry is given in Table El 

2.4. Comparison data 

The UVLGs in our sample, which have a median redshift 
of z ~ 0.15, lie at redshifts that are intermediate between 
local starburst galaxies (~ 0.02) and LBGs at high red- 
shift (z > 1.5). Throughout the paper, we therefore find 
it instructive to present our results with respect to the 
following two sets of comparison data: 

2.4.1. Archival imaging data of Haro 11 and VV 114 

As shown in H07 (see their Fig. 12), a number of very 
nearby blue compact starburst galaxies and (U)LIRGS 
fall very near to the boundary of the UVLG selection 
window as defined in the Lpuy vs. Ifuv plane. Because 
there is quite some leverage in the definition of a typical 
'LBG', we will compare our data to two of these objects, 
which have been found to possess some properties that 
are similar to high redshift LBGs. The two local compar- 
ison objects are Haro 11 and VV 114, both of which lie at 
z = . 02 (for details se e lKnop et al1ll994l: iScoville et all 
20001: iGliidader et all 12001 iBergvall fc Ostlinl [20021 : 
Kunth et alj 120031: iGrimes et alj l2006t l2007l and refer- 
ences therein). Haro 11 is a blue compact dwarf galaxy 
consisting of a number of UV bright knots believed to be 
in the process of merging. VV 114 is a merging system 
consisting of a UV-luminous Western component (VV 
114W), and an IR-luminous Eastern component (VV 
114E) with very little associated UV emission. 

In this paper, we make use of an ACS/HRC image of 
Haro 11 taken through the filter F220W (A c w2255A) ob- 
tained from the HST archive (Program 10575, PI: Goran 
Ostlin). For VV 114(W), we use a Space Telescope Imag- 
ing Spectrograph (STIS) image in the NUV (A C «2365A) 
from Program 8201 (PI: Gerhardt Meurer). 

2.4.2. High redshift samples from \Lotz et al\ H200& ) 

Our high redshift comparison samples consist of 55 star- 
burst galaxies at z ~ 1.5 selected from GOODS, and two 
(largely non-overlapping) samples of in total 82 z ~ 4 
LBGs in GOODS and the UDF from lLotz etall (f2006h . 
The galaxy sample at z ~ 1.5 is based on a large spec- 
troscopic sample of strong emission line galaxies. We 
used the publicly available GOODS and UDF data, and 
extracted postage stamps of the objects in -B435 for the 
z ~ 1.5 sample, and in Veo6+*775 for the z ~ 4 sample. 
A detailed desc ription of the dat a and sample selection 
can be found in lLotz et all (|2006[ ) . 
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3. RESULTS 

3.1. Images and notes on individual galaxies 

In Fig. [T] we show the GALEX FUV and NUV images. 
The UV, Ha, and optical images taken with HST are 
shown in Fig. [2] which have a resolution of ~ 50 times 
that of the GALEX images. Qualitative remarks on the 
morphologies of the UVLGs in Fig. [2] are as follows. 

SDSS J005527.46-002148.7. This object consists of a 
strong point source surrounded by diffuse emission that 
can be seen in C/330, zsso and Ha. The diffuse emission 
has a high surface brightness region extending to the NE 
with respect to the nucleus, which suggests either the in- 
fall of a small diffuse galaxy or an additional off-center 
star forming region. The continuum-subtracted Ha emis- 
sion is peculiar with an 'arm' of emission extending to 
about 1" South of the nucleus. The morphology of the 
Ha is very different from the zssa morphology, and at 
its outer extremity, the Ha is not associated with any 
UV emission. We suspect that the extended Ha is due 
to an outflow over the larger region probed by the ZS50 
continuum. 

SDSS J032845. 99+011150.8. The z S50 image sug- 
gests a merger of two low surface brightness galaxies, 
as evidenced by the tidal tails extending symmetrically 
from the center to the Northwest and Northeast. The 
starburst dominates the nuclear region, in which five iso- 
lated point sources can be seen. The easternmost knot 
dominates in Ha. 

SDSS J040208.86 050642.0. The z 850 image suggests 
a merger of two diffuse systems as the contours of faint 
^850 continuum emission show a sharp bend in position 
angle. The image center is dominated by three point 
sources, and several fainter ones lie further out on either 
side of the center. One of these faint, compact regions 
(about 1" South of the nucleus) is more extended in Ha. 
SDSS J080844.26+394852.4. There is a bright, un- 
resolved source seen in the ^330 image that is located in 
the centre of an extended, highly diffuse galaxy seen in 
Z85o. The unresolved source is remarkably bright as both 
U330 and Ha show the ACS diffraction spikes. The Ha 
image suffers from poor point source subtraction. The 
object seen ~ 2" to the South lies at a similar redshift, 
and may be interacting. The nucleus of this companion 
is also detected in L/330, but its flux is ~ 75 x lower. 
SDSS J092600.41+442736.1. The faint contours of 
the diffuse emission in z^o suggest that at least two dif- 
fuse galaxies may be merging. A bright, compact com- 
ponent dominates in ZS50 and Ha. The faint extension 
about 1" East of the nucleus in z 8 50 has little associated 
Ha emission. The 6^330 image failed to execute due to a 
guide star problem. 

SDSS J102613.97+484458.9. This galaxy has a ring- 
shaped morphology with most of the starburst occurring 
along the Eastern rim. There is also an isolated knot to 
the Southwest. The star- forming region contains several 
bright, unresolved components and diffuse emission. A 
very faint diffuse companion can be seen in zssoi about 
2" to the South, and very faint extended emission, pos- 
sibly tidal debris, is seen directly to the East along the 
entire extent of the galaxy. The overall morphology is 
very similar to that of the "drop-thr ough" ring galaxy 
NGC922 studied bv lWong et alJ (f2006h . 



SDSS J135355.90+664800.5. This UVLG is highly 
irregular in all bands, showing many star- forming knots 
and diffuse emission. The object is interacting or merging 
with a warped, edge-on or filamentary galaxy located just 
to the East, which is much redder and only just detected 
in C^330- A longslit spectrum along the major axis of 
the system shows that the two objects are likely to be a 
counter-rotating merger (Overzier et al., in prep.). 
SDSS J214500. 25+011157.6. This is the only object 
in our sample for which the ACS image provides a direct 
match at the wavelength of the GALEX FUV image. 
The F150LP image shows a very compact, but elongated 
starburst region in which we can discern three different 
knots. Fainter FUV emission comes from a slightly more 
extended region of The optical morphology in Veoe 
is strikingly different. The image is still dominated by 
the emission surrounding the starburst region, and a faint 
spiral structure may be present. A companion object 
appears to have gone straight through the star forming 
nucleus as evidenced by its tadpole-like morphology with 
a faint tail pointing towards the main galaxy and a trail 
of tidal debris that can be traced back to a location that 
is about equidistant from the nucleus on the opposite 
side of the main galaxy. 

3.2. Emission lines 

Although the sample of UVLGs does not contain any 
broad line AGN, we want to make sure that we are 
studying the morphologies of starburst galaxies rather 
than those of narrow line AGN. The main optical 
emission line diagnostic diagrams involving the lines 
[OIII]A5007A, H/3, [OI]A6300A, Ha, [NII]A6584A, and 
the [SII]AA6713,673lA doublet (measured from the SDSS 
spectra) are shown in Fig. [3] (see also H07). The 8 ob- 
jects observed by HST are indicated by the large filled 
circles and their IDs, while small filled circles indicate ob- 
jects in our supercompact UVLG sample without HST 
data. Most of the UVLGs lie along the main star forming 
sequence (points) in log([OIII]/H/3) vs. log([NII]/Ho) 
(left panel of Fig. [3]), albeit offset towards higher values 
of log([OIII]/H/3 probably due to a more intense ioniz- 
ing radiation field. It is interesting to note that a simi- 
lar offset from the main star-forming seq uence has been 
observed for high-reds hift LBGs as well (jShaplev et al.1 
[20051 : lErb et alJl2006aft . 

About one quarter of the sample falls i n the region 
between the relations of | Kauffmann et al.l (|2003f ) (solid 
line) and iKewley et al.l ()2006l ) (dashed line) that is typi- 
cally populated by objects having a composite spectrum 
consisting of a metal rich stellar population and an AGN. 
When we plot log([OI]/Ha) vs. log([OIII]/H/3) (middle 
panel) or log([SII]/Ha) vs. log([OIII]/H/3) (right panel), 
the objects that were in the AGN/starburst composite 
region move to the far, opposite side of the star form- 
ing track (see object SDSS J080844.26+394852.4 that is 
most relevant to this paper). It is not exactly clear what 
causes this behavior of the line ratios, but it is consis- 
tent with starburst mo del predictions hav ing very high 
ionization parameters (|Kewley et al.l |200TI ) as might be 
expected for our sample. A more detailed analysis of the 
ionization parameters and the possibly remaining contri- 
bution from AGN will be given elsewhere ( Overz ier et 
al., in prep.). In any case, Kauffmann et al. (2003) have 
shown that the continuum emission from local narrow- 
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line AGN is almost always dominated by its stellar pop- 
ulation; this assures us that the HST images accurately 
reflect the morphology of the stellar light. 

3.3. HST-based size measurements 

3.3.1. UV sizes 

Our sample of UVLGs is selected to have high surface 
brightness, which is estimated from the FUV flux and a 
half-light radius measured from fitting a seeing-convolved 
exponential profile to the SDSS u-band images (see Sec- 
tion 1, H05 fc H07). The resulting seeing-deconvolved 
sizes were on the order of 1-2 kpc, but were somewhat 
uncertain because most of the UVLGs are only barely re- 
solved in the w-band SDSS images. The F330W images 
are a factor 10 higher in resolution and thus provide a 
reliable and independent test for the true sizes of these 
UVLGs. We have used SExtractor's circular aperture 
method to estimate the radii; the 50 and 90% radii were 
computed using the light enclosed within 2.5 times the 
elliptical Kron radius also determined using SExtractor 
(Table [JJ . The sizes measured from the HST images are 
somewhat smaller than our previous estimates based on 
SDSS. In the case of SDSS J080844, approximately half 
the light in F330W is coming from a region that is still 
unresolved even at the resolution of the HRC (~ 0'.'07, 
or ~110 pc at z — 0.091). Excluding this extreme case, 
the radii range from 0.4 to 1.9 kpc. 

One could question (see Appendix) whether the sizes 
measured in C/330 or u-band (used to calculate the UV 
surface brightness) are representative for the sizes of 
the regions that emit in the FUV and NUV as seen by 
GALEX at the much lower resolution of ~ 4" (FWHM). 
We can test this simply by measuring the slope, /3, of 
the UV continuum (with f\ oc X 13 ) across the GALEX 
bands, and comparing the flux measured in the HST im- 
age to the extrapolated value at ^330. The UV slopes 11 
are given in Tabled and we note that they are very sim i- 
lar to that of high-redshift LBGs (e.g. lOuchi et alJl2007f ). 
The measured flux in C/330 is lower than the predicted val- 
ues by a factor of ~ 1.3 on average. The measured flux in 
u-band is ~ 1.2 x lower than the value extrapolated from 
the UV. The generally good correspondence between the 
predicted and measured continuum flux at 1500-3500A 
indicates that the emission structure seen in the HST 
and SDSS images is representative of the light distribu- 
tion in the FUV. A further test is provided by the case 
of SDSS J214500 for which we have an ACS/SBC image 
taken through F150LP, a filte r that is almost ide ntical 
to the GALEX FUV filter (see lTeplitz et alj |2006). The 
corresponding FUV surface brightness calculated from 
the ACS data is \og w Ipuv,Acs — 9.53 L© kpc -2 , in 
agreement with HO 7. 

3.3.2. Ha and Optical sizes 

The distribution of the optical and Ha light need not 
be the same. We find, however, that the sources are 
compact at all observed wavelengths (Table HJ), and the 
starbursts dominate the structures. The half-light radii 
measured from the continuum subtracted Ha images are 
slightly larger than those measured from [/330 (~0.5-2 
kpc). In some cases, the Ha morphology suggests an out- 
flow (SDSS J005527, and perhaps J032845 and J092600; 



see Fig. [2]), but this does not affect the measured sizes. 
The half-light radii in F850LP are typically 1-2 kpc, with 
the exception of SDSS J135355 which has r 50j 85o ~ 4 kpc 
(note, however, that it has a nearby companion). Both 
the 50 and 90% radii are typically twice as large as the 
UV size (~2-5 kpc), indicating that these systems have 
relatively faint, underlying structures that were already 
in place prior to the current episode of star formation. 

3.4. UV-optical Colors 

We have used SExtractor to derive (circular) t/330-2850 
radial color profiles out to the 90% flux radius. The 
results are shown in Fig. 2] Most objects have large 
gradients in their UV-optical color profile, with objects 
typically being bluer near the center, and redder further 
out. Because the profiles were calculated with respect to 
the centroid of the objects in the optical image, for some 
objects (032845, 040208 & 102613) the bluest regions 
are sometimes slightly offset from the nominal center of 
the galaxy (see Fig. [5]). We have defined an "inner" 
color measured within the t/330 half-lightradius and an 
"outer" color measured over the region between the 50 
and 90% radii in z§50 (see Table @|. In all cases the 
inner color is bluer than the outer color, and the steepest 
color gradients are seen for the most compact objects (in 
£^330) • This indicates that these objects are composites 
of recent, central starbursts within older (or more dusty) 
extended structures. 

3.5. Extinction 

For starburst galaxies there exists a good correlation be- 
tween the amount of UV flux that is absorbed by dust 
(usually taken to be a foreground "screen"), and the 



amount of flux (re-)emitted in the far-IR ( Meurer et al 



19971 ICalzettil [20011: IKong et alJ [200l ISeibert et al 
2005) . To estimate the internal extinction for our sample 
we have calculated the attenuation in the FUV, Apuy, 
using the ratio of the bolometric dust luminosity, Ltir, 
to bolome tric FUV luminosity, L puv, and the fitting 
formula of iBurgarella et all (|2005l ): 



Apuy = a\x + a 2 x + a 3 x + a 4 , 



(1) 



with [01,02,03,04] = [-0.028,0.392,1.094,0.546] and 
x = \o%w(L tirI 'Lfuv)- The bo lometric dust luminosity 
is given by (pale fc Heloul[2002]) : 



t>TIR = £l ^24-^24 + £2^70-^70 + £3^160^1 



(2) 



11 Calculated from Pgal 



■ 2.32(m FUV - m NUV )AB - 2. 



with [£1,6,6] = [1.559,0.7686,1.347]. Although we do 
not have observed flux densities at 160 /im, we can es- 
timate 160 /im using the strong empirical correlation 
between (hi hi) dust and (ho/fw o) by employing the 
synthetic models of I Dale fc Heloul (|2002t ). The corre- 
lation arises because more intense radiation fields have 
larger (hoi f wo) due to hotter large grains, and smaller 
(hi hi)dust due to larger emission by small grains at 24 
fim. The dust emission at 8.0 and 24 /jm was calculated 
from the observed flux densities using (/s)dust = h ~ 
0.232/ 3 .6 and (/ 2 4)dust = /24 - 0.032/ 3 .6, which includes 
a correction for the contributio n due to stars based on the 
3.6 fim flux pale et all 12005). Finally, we calculate the 
extinction using th e formula E(B — V) stais = Apuv I 'k(X) 
with k(X) « 10.77 (|Calzettill200lT ). The extinction ranges 
from E(B - V) « 0.01 to 0.14. The IR flux densities, 
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bolometric luminosities, Apuy and extinction are given 
in Table [5l The extinction will be needed in Section 4 
when we estimate masses and ages from the data using a 
se t of model star formatio n histories. We refer the reader 
to lBasu-Zvch et al.l (|2007t ) for a more detailed discussion 
and analysis of the extinction properties of this sample 
using alternative estimators. 

3.6. Morphologies 
3.6.1. Gini coefficient, M20, and Concentration 

In order to compare the morphologies of the UVLGs 
to those of nearby galaxies and of LBGs at high red- 
shift, we calculated the Gini coefficient (G; a measure 
of the equality with which the flux is distributed across 
a galaxy), M 2 o (the log of the ratio of the second or- 
der moment of the pixels containing the 20% brightest 
flux to the total second order moment), and concentra- 
tion (C; five times the log of the ratio of the circular 
radii containing 80 and 20% of the flux). To calculate 
the structural parame ters we follow the proc edures de- 
scribed in full detail in lLotz et all (|2004U2006f) . First, we 
use SExtractor to make an object segmentation map and 
mask out neighboring objects. The image is background 
subtracted, and we calculate an initial Petrosian radius 
(rp with r\ = 0.2) using the object center and (elliptical) 
shape information from SExtractor. We then smooth the 
image by a = rp/5 and create a new segmentation map 
by selecting those pixels that have a surface brightness 
higher than the mean surface brightness at the Petrosian 
radius. We recalculate the object center by minimizing 
the second order moment of the flux, and then recalculate 
the Petrosian radius in the original image using this cen- 
ter. The total flux is defined as the flux within a radius 
of 1.5 x rp and C is then calculated in circular apertures. 
The individual measurements are listed in Table [6l and 
the results are shown in Fig. [5l For each object, the 
C/330 measurements are plotted as stars and the zsso as 
circles, and are connected by a dotted line. For 080844, 
G and C values in C/330 must be considered lower limits, 
and its M20 value is an upper limit because the central 
pixel of the unresolved component contains «20% of the 
total flux. 

Although our sample is small, we will attempt to 
make a general comparison with the morphologies of 
nearby galaxies. As shown in Fig. [5j the UVLGs pop- 
ulate the region roughly defined by M 2 q < —1.3 and 
0-53 < G < 0.7 (left panel), and 3 < C < 5 (right 
panel). We can compare these measurements to the vari- 
ous diyisionsof morp hological parameter space discussed 
by lLotz et al.l (|2006h . Our UVLGs have smaller M 20 
than merging galaxies with two clearly separated nuclei 
(M20 > —1-1, left hatched region), larger G than the ex- 
pectations for exponential disks (dotted box), M 2 o sim- 
ilar or larger than those expected for bulge-dominated 
objects (right hatched region), and G and C that are 
either similar or smaller than those expected for bulge- 
dominated galaxies seen face-on (solid box). We do not 
find any systematic differences between the morphologies 
measured in C/330 and Z850 for the sample as a whole. 
The morphological parameters in the UV and optical 
are thus dominated by the same (bright) regions, even 
though some of the objects possess very faint, extended 
structures in zsso- 



Fig. [5] shows that two UVLGs qualify as bulge- 
dominated in C/330, while four objects qualify as bulge- 
dominated in zg5o. None of the objects qualifies as a 
pronounced merger, and only one object lies just on the 
border of the exponential disk region in zsso, with the 
remainder populating the regions in between. The over- 
all distribution of morphologies is very similar to that 
of star forming galaxies at z ~1.5 an d LBGs at z ~ 4 
selected from GOODS and the UDF. lLotz et alJ (|2006l ) 
find median values of G ~ 0.55, M20 ~ —1.5 and C ~ 3.3 

atz- 1.5, and G ~ 0.58, M 20 1.6 and C ~ 3.8 at 

z ~ 4 (encircled crosses in Fig. |5|). However, a fair com- 
parison of the morphologies of UVLGs with the higher 
rcdshift comparison samples should be carried out at the 
same S/N and resolution. In the next section, we will 
carry out such a comparison. 

3.6.2. Redshift simulations 

We now investigate whether the conclusions of the pre- 
vious subsection are still upheld if the analysis is carried 
out using images that are simulated to have the same 
depth and resolution as LBGs observed at high redshift. 

We fo llow common prac t ice (e . g.. see | Giavalisco et alJ 
I1996U Ifflbbard k. Vaccal 11997k TBouwens et ajj 119981 : 
iPapovich et "all 120031 : IConselicd l2003t lLotz et al.ll200l . 
and apply corrections for cosmological surface brightness 
dimming and for changes in physical resolutions to the 
C/330 images, in order to match the instrumental con- 
ditions under which our objects would be observed in 
typical surveys when placed at z = 1 — 4. 

The first step of the procedure is to rebin the C/330 
images by a factor b = {0\/ 9 2 ){s 2 / s\), where 0j is the 
angle on the sky of an object of fixed size d at z = Zi, and 
Si is the instrumental pixel scale (in arcsec pixel -1 ). The 
rcbinning factor can be expressed in terms of either the 
angular diameter distance, D A t = d/9i, or the luminosity 
distance, Dp i — (l + Zj) 2 !)^, at low and at high redshift 
(z 2 > 21): 

, = Da 2 S2 = ( 1 + zi \ 2 D L2 s 2 ( ) 

D Al si \l + z 2 j D Ll V 1 ' 

The second step is to reduce the surface brightness of 
each (rebinned) pixel according to the relative amount 
of cosmological dimming of a galaxy at z 2 with respect 
to that at z\. We calculate the scaling by making use 
of the fact that the absolute rest-frame magnitude (or 
luminosity) of the object before and after redshifting will 
be conserved (/Wa 2 /(i+2 2 ) = -^/ai/(i+zi) with matched 
filters so that \ 2 = Ai(l + z 2 )/(l + zi)). 

In order to compar e our results to the high redshift 
samples of lLotz et al.l (|2006f) . objects were simulated at 
z = 1.5 in -B435, and at z = 3.0 and 4.0 in both Vqq§ and 
2775 . Objects were simulated at the depths of GOODS (3, 
2.5 and 2.5 orbits in -B435, V606 and 1775, resp.) and the 
UDF (56, 56, 150 orbits in -B435, Vqq% and 1775, resp.). For 
completeness and future reference, we also simulated the 
sample at the shallower COSMOS survey in the isi4 filter 
(1 orbit). For each filter and for each redshift, we used 
the HST images that correspond most closely in terms 
of rest-frame central wavelength. Where necessary, we 
rescaled the flux of the input images assuming the UV 
slopes determined in Section 13.51 The final steps of the 
simulation consist of applying Poissonian noise to the 
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simulated profiles based on the typical exposure times 
of COSMOS/GOODS/UDF, convolving the images with 
a Gaussian to match the desired output PSF size (~ 
C//12), and placing the simulated object inside an empty 
region in a COSMOS, GOODS or UDF image to obtain a 
realistic background. The simulated pixel scale was 0'.'03 
pixel" 1 for GOODS and the UDF, and 0"05 pixel" 1 for 
COSMOS. 

The rest-frame UV images of Haro 11 and VV 114 
(see top panels of Fig. O were artificially redshifted in 
an identical manner using a two-step process. First, we 
simulated their F330W images at a redshift (z = 0.15), 
depth (2500 s), plate scale (07025 pixel" 1 ), and seeing 
(~ 0'.'075) comparable to our observations of the UVLGs. 
These images are shown in the bottom panels of Fig. [6l 
Both Haro 11 and VV 114 possess multiple bright nuclei 
apart) when observed in U330 at z = 0.15. VV 
114 also shows a significant amount of diffuse emission 
in between the knots. Next, these images were used to 
simulate how these objects would appear at high redshift, 
analogous to the UVLG simulations described above. To 
match the desired rest-frame wavelengths of the output 
filters, small flux extrapolations were performed using 
= — 1-4 f or both objects a s mea sured bv lGoldader et al.l 
(12001 and lBergvalleTall (f2006h . 

3.6.3. Results 

Figs. [7J [51 and [H compare postage stamps of the local 
LBG analogs (including Haro 11 and VV 114) that are 
simulated at the depths of the GOODS, UDF, and COS- 
MOS surveys. Results are shown at the three different 
redshifts discussed above (z = 1.5, 3.0 and 4.0). The im- 
ages demonstrate that as redshift increases, components 
that are well separated at low redshift blend due to the 
lower spatial resolution, and that low surface brightness 
features are lost due to surface brightness dimming. The 
latter effect is most severe in the shallow COSMOS data, 
and least severe in the deep UDF data. 

Fig. [TD] shows how the measured morphological pa- 
rameters G (top panel), M20 (middle panel) and C (bot- 
tom panel) change as a function of redshift and survey 
depth. There is a systematic drop of ~0. 05-0. 10 in G 
from z ~ 0.15 to z ~ 1.5, followed by a further less sig- 
nificant decrease of a few hundredths out to z = 4. The 
first decrease is caused by the strong drop in resolution, 
as well as the loss of faint extended features that cause 
the flux to be more ev enly distributed o ver the galaxy 
profile. Simulations bv lLotz et a l. (2006) show that the 
second decrease is mainly caused by the loss of low sur- 
face brightness pixels lowering G somewhat further at 
z > 1.5. The loss of low surface brightness features with 
redshift also tends to increase M20 (the total second or- 
der moment is lowered while the second order moment 
of the 20% brightest flux stays roughly constant), and 
tends to lower the concentration, because the 80% flux 
radii are systematically underestimated. The measure- 
ments at the COSMOS depth show quite a large scatter 
with respect to GOODS and the UDF, indicating that 
the lower S/N and its larger pixel scale result in rela- 
tively unstable morphology measurements for LBGs. 

The main results of our morphological comparison are 
presented in Fig. [Til We compare the UVLG morpholo- 
gies (filled stars) with those of the low redshift LBG 
analogs (open stars) and the high redshift comparison 



sample (plusses) of lLotz et al.l (|2006). In order to make 
sure that the morphologies are measured in a consistent 
manner, we used our own code to recalculate the mor- 
phologies of the 55 z ~ 1.5 starburst galaxies in GOODS 
(bottom panels), and the 82 z ~ 4 LBGs in GOODS 
(middle panels) a nd the UDF ( top pa nels) analyzed by 
lLotz et all ()2006[ ). lLotz et al.l (|2006f) showed that the 
morphologies of the z ~ 1.5 sample are very similar to 
that of the z ~ 4 sample. The distribution of the UVLGs 
is very similar to that of the high redshift objects as well, 
although the latter has a larger scatter due to the much 
larger sample size. 

3.6.4. Comparison with Haro 11 and VV 114 

Haro 11 and VV 114 (open stars in Fig. [TTj) have higher 
M 2 o and lower G and C when compared to the UVLGs. 
We can understand this given that both objects possess 
several bright nuclei with a separation that enables them 
to be significantly resolved even when simulated at z = 
0.15 — 4.0. Their qualitative morphology perhaps most 
closely resembles t hat of the high redshift "clump-cluster 
galaxies" fe.g.. see lElmegreen et al.l l2004'). Although ob- 
jects with similar morphological characteristics are cer- 
tainly present in high redshift samples, they do not make 
up the majority of LBGs as can be seen from the distri- 
butio n in morphological types of LBGs (e.g. lLotz et al.l 
120061 iRavindranath et all [2001 lElmegreen et all 120071 . 
and Fig. [TTJ in this paper). 

It is important to note that our GALEX/SDSS selec- 
tion of local LBG analogs at z « 0.1 — 0.3 is likely biased 
against finding objects having such widely separated nu- 
clei as Haro 11 and VV 114. Because we selected objects 
principally on having a high FUV luminosity, as well as 
a high UV surface brightness based on the half-light ra- 
dius measured in the w-band, we are most sensitive to ob- 
jects that are not or barely resolved in SDSS. The surface 
brightness requirement is necessary in order to obtain a 
relatively clean separation between the LBG analogs on 
one hand, and large, UV-luminous (predominantly spi- 
ral) galaxies on the other. However, we will further in- 
vestigate whether we can find more LBG analogs of the 
type of Haro 11/ VV 114 by looking closer at the objects 
that straddle the boundary of the LBG analog selection 
criteria in a future paper (Overzier et al., in prep.). 

3.6.5. Summary 

We conclude that the morphologies of the UVLGs can- 
not be distinguished from the morphologies of high red- 
shift LBGs when measured from the redshifted images 
simulated at the same depth as GOODS and the UDF. 
This further strengthens our conclusion that the rela- 
tively nearby supercompact UVLGs and high redshift 
LBGs are very similar. 

In the following sections, we will turn our attention 
again to the undegraded HST images of Fig. [2l In Sect. 
[4] we will study in detail the nature of the starburst re- 
gions, and in Sect. Owe discuss the connection between 
morphology and the mechanisms responsible for trigger- 
ing the vigorous star formation observed in these local 
LBG analogs. 

4. SUPER STARBURST REGIONS 

As discussed in Section 3, in all our sources the UV light 
is distributed in a series of bright, unresolved knots em- 
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bedded in a region of more diffuse emission. In this sec- 
tion we will show that these regions can be interpreted as 
being (super) starburst regions (SSBs). An understand- 
ing of the physical nature of these regions is crucial for 
interpreting the UV emission from LBGs. 

4.1. Identification and photometry 

In Fig. [12] we indicate all the starburst regions iden- 
tifiable by eye from the E/330-band images (small cir- 
cles). Although our identification by eye is somewhat 
subjective, most regions are bright and isolated and thus 
easily recognised. We identify a total of 41 regions in 
our sample of 8 galaxies. Some galaxies contain only 
one knot, while others contain as many as ten. Next, 
we measured the fluxes and colors of each region us- 
ing circular aperture photometry in matched C/330 and 
z 850 images (resolution of ~ 0'.'12, FWHM). We were 
able to isolate all of the starburst regions identified in 
Fig. [T2] using circular apertures of Of! 3 in diameter (cor- 
responding to physical radii of ^240- 480 pc). The encir- 
cled energy (EE) measurements of (jSirianni et al.ll2005f ) 
show that circular apertures of this diameter enclose 
~ 80% of the total light of a point source observed with 
ACS. We did not apply the EE correction. Magnitudes 
were corrected for Galactic extinction, and the magni- 
tude distribution of the 41 starburst regions is shown 
in Fig. [13] Absolute magnitudes were calculated from 
M = m - 51og 10 (D L (pc)) + 5 + 2.51og 10 (l + z). We 
did not apply a if-correction as the central rest-frame 
wavelengths at which the absolute magnitudes are deter- 
mined differ only by a small amount (< 300 A). Fig. [13] 
shows that our 'by eye' selection is relatively complete 
at TO330 < 23 mag (-^330/(1+*) ^ —16.5). We also deter- 
mined the flux and color of the large, diffuse region that 
generally surrounds the compact starburst knots (large 
circles in Fig. [T2]) . These regions measure ~2-5 kpc in 
radius. The fluxes measured in the small apertures were 
subtracted from the larger aperture. 

We find that the UV light is typically dominated by 
emission from the compact, luminous regions: the com- 
bined flux of the unresolved regions contributes an aver- 
age of ~50% to the total flux in [/330. Object 102613 has 
the lowest contribution from SSBs (~ 34%) and object 
080844 the highest (~ 80%). Although we do not have 
a t/330 image for 092600, we can infer that its structure 
is very similar to the other objects given that the unre- 
solved knots contribute ^42% to the total flux in zsso, 
and its morphology in Ha (assumed to be a good proxy 
for its C/330 morphology, see Fig. [2]), shows a similar 
knotty structure. 

4.2. Color-magnitude diagram 

In Fig. rn] we show the color-magnitude diagram of all 
regions identified in Fig. [12] Objects 092600 and 214500 
were omitted. Starburst regions from each galaxy are 
plotted using a different filled symbol, with their corre- 
sponding ID number from Fig. [I2]indicated for reference. 
Open symbols of corresponding shape indicate the diffuse 
regions that surround the compact regions. In order to 
compare the measured colors and magnitudes to model 
starburst tracks, we 'dereddened' the measurements us- 
ing the global reddening value for the stellar continuum 
that was derived from the bolometric dust to FUV lumi- 
nosity ratio (see Sect. 13.51) . We note that this correction 



is in general quite small. The compact regions span quite 
a large range both in color and in absolute magnitude 
(filled symbols in Fig. [Ll| . The diffuse regions, defined 
by the larger areas in Fig. [12] but with the circular re- 
gions around all the identified point sources removed, lie 
at Mjj ~ —20 and {Us^q-zs^o)^ 0.5 (open symbols in 

Fig. E]). 

To inter pret these colors, w e have used STAR- 
BUR ST99 (jLeitherer fe Heckmanl Il995t iLeitherer et~aTl 
1999) to predict the color and magnitude evolution of 
a starburst as observed through our filters. The tracks 
in the left panel of Fig. [14] represent an instantaneous 
burst model with burst masses of M = 10 7 Af© and 
M = 10 s Mq, while tracks in the right panel indicate 
a continuous star formation model having SFRs of 0.1, 1 
and 5 M© yr _1 . We use a Kroupa initial mass function 
(IMF) with a slope of a = 1.3 between 0.1 and 0.5 M© 
and a = 2.3 between 0.5 and 100 M©. We also use the 
Padova 1994 models with thermally pulsating AGB stars 
added (I Vazquez fc Leithererl [20051 1. The tracks shown 
were obtained after redshifting the rest-frame spectra to 
z = 0.15, the average redshift of our sample, and mea- 
suring the magnitudes in U330 and zsso- The blue lines 
correspond to a metallicity equal to that of the Large 
Magellanic Cloud (Z = 0.008), while the red lines are 
for solar metallicity. Ages (in Myr) have been indicated 
along the tracks. Dashed lines indicate the contribution 
from the nebular continuum (assuming no reddening). 
The general behaviour of all tracks is to become redder 
with age. After about 10 Myr, the colors become steadily 
redder, because the most massive O stars have become 
supernovae and lower mass O stars have become red su- 
pergiants. The red supergiant dominated phase at 10-20 
Myr is highly dependent on metallicity, as their number 
is lower and their temperature is warmer in l ower metal- 
licity starbursts (] Vazquez fc Leithererl [20051 ). After 20 
Myr, O stars have disappeared leading to a gradual red- 
dening with age. Continuous models become brighter, 
while instantaneous models fade with time. 

4.3. Results 

The instantaneous burst and continuous star formation 
models have been chosen to bracket the range of pos- 
sible star formation histories of our compact starburst 
regions. Most of the regions are clearly detected in Ha 
indicating that they must at least have had very recent 
star formation. In Fig. [15] we plot the SSBs in the plane 
of color vs. the rest-frame equivalent width (EW) of Ha 
measured for each SSB. The measured equivalent widths 
are 30-500A, and are seen to correlate with the (7330-2:850 
color, suggesting that the Ha EW is a reasonably good 
age indicator. We have overplotted the predicted EW of 
Ha based on the STARBURST99 models shown in Fig 
[T4l For instantaneous models, the Ha EW is less than 
~100A after 6 Myr, and less than 10A after only 10 Myr 
(dashed lines). For continuous star forming models, a 
high Ha EW is maintained over a much longer period of 
time, with EW Ha >100A for ages <1 Gyr. 

Comparison of the instantaneous burst and continuous 
star formation tracks with the observed compact and dif- 
fuse regions lead us to the following conclusions: 
(i) The M ~ 10 7 ~ 8 M© instantaneous starburst tracks 
shown in the left panel of Fig. [H] provide a generally 
good match to the majority of the compact regions, but 
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the ages of ~10-100 Myr inferred are inconsistent with 
their large Ha EWs, as shown in Fig. [T5j A better match 
to both the color-magnitude diagram and the Ha EWs is 
provided by the continuous star formation models shown 
in the right panel of Fig. [TU The exact ages are uncer- 
tain, because of small-scale differences between the C/330 
and Ha morphologies (note that the apertures were de- 
fined on the C/330 image), uncertainties in the small-scale 
dust distribution (now assumed to be a single foreground 
screen), and because of the unknown fraction of the light 
that is due to the background population inside each of 
the small apertures (not subtracted). Our best estimates 
for the ages of the starburst regions range from a few tens 
of Myr to a few hundreds of Myr based on a comparison 
with the tracks in Fig[l4j These ages should be regarded 
as upper limits, because we may have underestimated 
the contribution from dust and older stars. The scatter 
in Mjj magnitude can be interpreted as a sequence in 
SFR ranging from ~0.1 to ^5 M© yr _1 (or total burst 
mass ranging from ~ 10 7 Mq to several 10 8 Mq). The 
scatter in (J/330— £85o) color can be interpreted as a se- 
quence in age (or burst age in the case of instantaneous 
bursts). Note that our results are almost independent of 
metallicity. 

(ii) For the brightest starburst regions, there is a strong 
tendency for neighbouring regions to have a very simi- 
lar color and magnitude, e.g. regions (2,4,7) in object 
135355, regions (1,4,5) in object 040208, and regions 
(4,5) in object 032845. This suggests that the starburst 
in these objects tends to occur in the form of several 
co-eval and equally massive bursts. 
(Hi) The two compact regions that are outliers (in the 
sense that they appear to be both extremely young, ~6 
Myr, as well as very massive, M* ~ 10 8 Mq) each corre- 
spond to objects where the UV emission is dominated by 
a single, luminous knot (objects 005527 & 080844). Both 
objects have a very large Ha EW, indicating that they 
are massive, young starbursts. The SDSS spectrum 12 of 
005527 shows a broad feature around a bright and nar- 
row Hell 4686A line that is characteristic of Wolf-Rayet 
(WR ) galaxies (e.g. IContil 119911 : iSargent fc Filippenko 
119911: llzotov et al.lll997D . The WR phase is a relatively 
short phase in the life of the most massive stars during 
which they experience a large mass loss. The presence 
of WR stars is a unique indicator of young starbursts, 
as WR stars disappear after ~ 6 — 10 Myr depending on 
the metallicity of the burst, with less mass loss oc curring 
at lower metallicity (jLeitherer &: Heckmanl[l9 95). Fur- 
ther confirmation of this young age comes from the large 
equivalent width of H/3 (~60A), consis tent with an age 
of -6 Myr (jLeitherer fc Heckmanlll995h . 

(iv) Some other regions with very blue colors and high 
EW in Ha (e.g. region 2 in 032845 and region 1 in 
102613) appear to be as young as the previous two, but 
are ~2-3 magnitudes fainter. 

(v) The more diffuse, inner regions in which the SSBs are 
embedded (open symbols in Fig. I14p are typically a few 
tenths of a magnitude redder than the average color of 
the compact SSBs they surround, and their total mag- 
nitudes are very similar to the combined flux of their 
SSBs. These regions contain an older stellar population. 
Given their large physical extent, the instantaneous burst 



model is unlikely to be appropriate. We find that the col- 
ors and luminosities of these regions are in good agree- 
ment with continuous star formation models with SFRs 
of ~ 1 - 5 Mq yr" 1 and ages of ~ 0.2 - 1.0 Gyr. 
(vi) The outer annuli (defined as the region between 
r 50,850 and r 90,850i see Tabled]) are consistent with older 
populations having ages between a few Gyr to a Hubble 
time, depending on the star formation history. 

4.4. Summary 

The starburst regions observed in the local LBG analogs 
are highly compact (radii of 100-300 pc), and charac- 
terized by bright, unresolved knots of emission within 
a larger region of diffuse star formation that extends 
out to a radius of a few kpc. The total UV emission 
from these compact starburst regions is substantial, in- 
dicating that 30-80% of the total SFR is generated in 
these knots. Comparison with STARBURST99 evolu- 
tionary tracks indicates that some of the regions are due 
to very young (~6 Myr) bursts, while other regions may 
have been forming stars in a more continuous manner for 
several tens to several hundreds of Myr as indicated by 
their relatively red colors and large Ha equivalent widths. 
There is a tendency for neighboring knots to have simi- 
lar colors and luminosities. This implies that they may 
be co-eval and of similar mass. Most of the SSBs are 
still unresolved in the unbinned U330 images (« 0'.'075, 
FWHM). The stellar mass densities implied are ~ 10 2-3 
Mq pc -2 estimated from their masses of ~ 10 7-8 M© 
and assuming an effective radius of ^100 pc. 

The masses of the super starburst regions are one to 
two orders of magnitude larger than the masses of the 
most massive clusters found in the local Universe. We 
may therefore conclude that the (unresolved) super star- 
burst regions s een in Fig. [T2| are likely composed of 
smaller units. iMeurer et al.l (|1995l ) found that typical 
starburst regions in local starburst galaxies consist of dif- 
fuse stars as well as clusters in the ratio of respectively 
80 to 20%. The most prominent of the clusters, the so- 
called super star clusters (SSCs), have masses of w 10 5-6 
Mq and sizes of a few to ~10 pc and likely correspond 
to globular clusters in the process of formation (provided 
the stellar IMF extends down to 0.1 Mq). 

We can get an even better idea of the sub-resolution 
structure of the super starburst regions by comparing 
the starburst regions in Haro 11 and VV 114 shown in 
the top panels of Fig. [6] with the same regions simu- 
lated at z = 0.15 shown in the bottom panels of the 
same figure. The majority of the bright nuclei seen in 
the images simulated at z = 0.15 are a blend of numer- 
ous smaller super starclusters that can clearly be seen in 
their unredshifted images. However, the most Northern 
and brightest knot in the redshifted image of Haro 1 1 is 
still largely unresolved in its unredshifted image. This in- 
dicates that single, highly massive star clusters are capa- 
ble of dominating the rest-frame UV/optical morphology 
of LBG-like galaxies at virtually any redshift. A detailed 
determination of the sizes, masses and ages of the star 
clusters in Haro 11 and VV 114 is beyon d the scope of 
this paper, and c an be found elsewhere (jScoville et alJ 
2000: iHayes et alJl2007t Adamo et al, in prep.). 

5. DISCUSSION 



12 http://cas. sdss.org/dr5/en/tools/explore/obj. asp?id=588015508738539596 5 ' 1 ' Triggering mechanism 
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Although the number of objects in our HST sample is 
currently quite small, the results presented in this paper 
can nevertheless guide us in the interpretation of the 
morphologies of these local starburst galaxies, and by 
extension, those of LBGs at higher redshifts. 

(i) Are the morphologies evidence for merging? 



Although some of the SDSS images showed evidence for 
close companions or faint extended emission, the galaxies 
were extremely compact and it was not at all clear from 
the SDSS images that they were highly disturbed. The 
biggest surprise from the HST images was that most of 
the UVLGs show disturbed morphologies at scales well 
below the SDSS seeing or sensitivity. Each UVLG shows 
at least one of the typical signs of merging or interac- 
tion, such as multiplicity or twists in the faint, outer 
isophotes (e.g., 005527, 032845, 040208, 092600), tidal 
debris or tails (e.g., 032845, 102613, 214500), and close 
companions (e.g., 080844, 102613, 135355, 214500). 

While none of the UVLGs appear as fully-formed Hub- 
ble sequence galaxies, in one case, 214500, a small spi- 
ral structure seems to be present both in the FUV and 
in V(306- Nonetheless, even in this object there is cir- 
cumstantial evidence that the starburst was triggered 
by a merging event, as evidenced by its companion that 
appears to have gone straight through its nucleus leav- 
ing behind a trail of tidal debris. The projected dis- 
tance between the centroid of the starburst region and 
the companion is about 5 kpc. The companion galaxy 
would traverse this distance in ~10-100 Myr assuming 
a velocity of a few hundred km s _1 and a straight tra- 
jectory in the plane of the sky. If we asume that the 
starburst is triggered by the interaction (as suggested 
by the VI506 morphology), this timescale constrains the 
maximum possible age of the burst. The results are 
consistent with 30-100 Myr ages derived from the col- 
ors and magnitudes of these systems. In many of the 
other cases, the region that dominates the UV morphol- 
ogy often lies at the interface of two diffuse merging 
structures, as identified by the contours of their outer 
isophotes in zsso- In the remaining cases, a compan- 
ion is always seen at a distance of less than <5 kpc 
from the UVLG. It is well-known that close neighbours 
can affect the SFR in galaxies, with close pairs hav- 
ing separati ons of < 50 kpc showing enhanced star for- 
mation (e.g.lLarson fc T inslev 1978; Lambas ct al. 120031 : 
iNikolic et aI]l2004HWoods et al.ll2006l : IOwers et~aHl2007l 
Li et al. in prep.) This supports the idea that the star- 
bursts are linked to mergers. The merging galaxies must 
be very gas-rich in order to trigger the very high level 
of enhanced star formation and the luminous super star- 
burst regions that we observe. 

It has been proposed that starburst galaxies at 
high redshift, such as Lyman break galaxies, are trig- 
;ered by mergers of relatively gas-rich disk galaxies 
Somerville et al.ll2001f ) . Simulations indicate that merg- 
ers of gas-rich galaxies can trigger massive starbursts. 
The conditions of the starburst depends on a num- 
ber of key parameters of th e merging system (e.i 



to the center is generally not disrupted or depleted due to 
star formation during the early stages of the merger, al- 
though the inflow of gas can be limited due to resonances 
if a bulge component is present, reducing the central gas 
densities and overall SFR and duration of the central 
starburst. The presence of gas is not always a sufficient 
condition for generating a starburst. The star formation 
efficiency of a merger depends largely on the galaxy spin 
direction (retrogade encounters have higher efficiency), 
and on the amount of gas that is be ing expelled due to the 
tidal interaction at first pa ssage (|Springel fc Hernquistj 
[20051 Idi Matteo et alJl2007h . 

We remind the reader that typical (L~L* =3 ) LBGs 
are not dwarf galaxies, but have masses of ~10 10 M© 
and SFRs that are much higher than those of t ypical 
(unobscured) local starbursts (jPapovich et al.l 120011 : 
Shapl ev et al.1 12001|). They do not s how evidence for 
spiral structure (|Papovich et al.ll2005f) . As such, they 



are very similar to our local LBG analogs. We conclude 
that the main mechanism responsible for triggering star 
formation in both high redshift LBGs and their local 
analogs is likely to be mergers. 

(ii) Are the compact super- starburst regions triggered by 
the mergers? 

Luminous star clusters (e.g. SSCs) occur in many ir- 
regular and dwarf galaxies (e.g. 30 Dor in the LMC). 
The triggering of compact SSBs as massive as the ones 
in our sample of UVLGs probably requires a highly 
efficient inflow and compression of the gas that is 
typically only seen in galaxy mergers or interactions. 
The inflow ensures that the pressure in the interstel- 
lar medium becomes larger than the internal pressure 
of giant molecular clouds, and that the clouds can col- 
lapse and form SSC s before the y are disrupted by su- 
pernova explosions (jBekkil 120041 ). SPH simulations of 
galaxy mergers suggest that the properties of the star- 
burst are determined by the magnitude, timescale and 
geometry of this inflow (se e Mihos & Hernauist 1994al : 
Bekki fc Couchll2001tlBekkill2004Hdi Matteo et al.ll2007h . 
Mihos fc Hernauist! (|1994af ) found that the central star- 



Mihos fc Hernauist 1994a. a; ISpringel fc Hernquist 



.g. see 
-d l2005b 



di Matteo et alJl2007L and references therein). During a 



merger of disk galaxies, bar-like patterns develop that ef- 
ficiently funnel gas into the central region. The gas flow 



bursts resulting from gas-rich mergers can have total 
sizes of ~ 350 pc, and contain most of the total mass 
in new stars. iBekki fc Couch! (|2001l ) predict that (i) 
SSCs have a narrow age distribution (~ 200 Myr) be- 
cause they form most efficiently during the stage when 
the gas density is the highest, (ii) the total number and 
intrinsic masses of SSCs are larger in major mergers due 
to larger tidal effects and higher peak SFRs, and (hi) 
SSC prod uction is more efficient in merging galaxies (e.g. 
ULIRGS. [Sanders fc Mirabe]||1996D th an in tidally inter- 
acting galaxies (e.g. M82. lKennicuttj|1998l ). 

Observations show that starbursts indeed occur in re- 
gions where the gas de nsities are the high est, typically in 
the cores (e.g. Arp 22Q. lShava et al.|[l994h and in circum- 
nucle ar starburst regions (e.g. VV114, iScoville et all 
I2000T ) . Stars may form in large complexes having to- 
tal stellar masses of ~ 10 7-9 M@, and containing sev- 
eral tens to hundreds of SSCs as well as smaller star 
clusters. This is true both in ULIRGS and in compact 
blue galaxies fe g.lMeurer et alj|1995t lOstlin et alj|1998t 
IWhitmore et al.l 119991 : IZepf et alJ Il999t iGoldader et all 
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120021 and therein). 

We conclude that the characteristic morphologies, 
ages, sizes and masses of the starburst regions we have 
observed are consistent with those predicted by simu- 
lations of gas-rich mergers. They are also analogous to 
the large star-forming complexes observed in very nearby 
merging galaxies. 

5.2. Low redshift lessons for high redshift LBGs 

Because of their relative proximity, studies of local star- 
forming galaxies can help us understand the nature of 
starburst galaxies at high redshifts. As explained in Sec- 
tion 1, up to now these comparisons have been of limited 
usefulness, because typical local starburst galaxies are 
systematically different from LBGs at high redshift. The 
sample of local LBG analogs first studied by H05 was 
specifically designed to bridge this gap. In this paper, 
we have studied the detailed morphologies of these sys- 
tems and we now discuss how our results impact on the 
study of star- forming galaxies at high redshifts. 

5.2.1. Morphologies and sizes 

Our analysis of the HST data has established that 
UVLGs have a light distribution dominated by very com- 
pact starbursting regions (ranging in number from 1 to 
more than 10 per galaxy) that are triggered by mergers 
or interactions. In most cases, an older, diffuse, stel- 
lar component is also present. Our redshift simulations 
showed that the sizes and morp hologies at_ J g_ = 1.5 — 4 



are similar to that of LBGs (e.g. iGiavalisco et alj 



Conselicd 120061: lLotz et all 120061: lOverzier et"aH 120071: 
Younger et al. [" 120071) . The morphological parameters 



1996a! 



measured at rest-frame UV and optical wavelengths are 
very similar. This is unlike typical local galaxies, but 
very simil ar to LBGs as both are dominated by young 
stars fe.g. iPapovich et aTll2005D . 

The redshift simulations presented in Section 3 yield 
a number of important conclusions: (i) at high redshift, 
the compact starburst regions can blend giving rise to a 
galaxy image that appears to have the one or two bright 
'knots' with a total half-light radius of 1-2 kpc, (ii) the 
diffuse emission is partly lost in the background noise 
or may blend with the knots, and (iii) the undegraded 
optical images were in most cases required to establish 
a definite connection between the (UV) starbursts and 
merging (Section 5.1). 

If 'knots' are present in LBGs they are typically 
not resolved in HST observations, ground-based spec- 
troscopy _with__ada^ive__optfos_j(resolution ~ O'.'l or 
~1 kpc: iGenzel et alll2006t lLaw et aLll2007b[ ). or stan- 
dard spectroscopy (resolution ~ 0^5 or —4 kpc; 
IForster Schreiber et al]|2006t lErb eT al. 2006b). Because 
imaging observations of high redshift LBGs lack the 
sensitivity and resolution to see the (often subtle) fea- 
tures associated with merging, 2D or 3D spectroscopy 
is the most promising method of studying the connec- 
tion between star formation and morphology in these 
objects. Small samples studied to date show complex 
kinematics and high nuclear gas f ractions that are consis- 
tent with the merger hypothe sis (lErb et al.ll2003l I2006U 
IForster Schreiber et all 120061 : lLaw et alj l2007bQ . How- 
ever, some star-forming objects at z ~ 2 may have more 
order ed underlying structures (jForster Schreiber et al.l 
2006; Gcn zel et all 2006: W right et al.ll2007f) . 



Coming back to the important, long-standing question 
of whether the irregular morphologies of LBGs are a sign 
of merging or of pa tchy star formati on within a single 
(forming) disk (e.g. lLaw et all l2007a[ ) . it is safe to say 
that in virtually every single case presented in this paper 
we do not see any direct evidence for merging based on 
the HST UV images alone. However, in virtually every 
single case evidence for merging is readily apparent from 
the rest-frame optical image, and we reiterate that most 
of the features suggestive of merging would be too faint or 
too small to be seen at z— 2-4. In some of the local LBG 
analogs the burst occurs predominantly in only one of the 
members of a merging or interacting pair. In others the 
merger seems in such an advanced state that one cannot 
speak about separate systems any more, making the dis- 
tinction between 'merging' and 'patchy disk' essentially 
irrelevant. We conclude that the UV morphologies of 
LBGs are likely patchy as the result of merging. 

Although both the UV and optical morphologies are 
dominated by compact starbursts, some of our ob- 
jects show evidence for an older, extended component, 
which is generally no t seen in high redshift LBGs (e.g. 
IPapovich et al.ll2005h . In this respect, the UVLGs may 
be more evolved than LBGs at high redshift. However, it 
is important to remember that the age of the universe at 
z ~ 3 is only ~ 2 Gyr. The spectral energy distributions 
of LBGs are consistent with the presence of such a "max - 
imally old" stellar component IPapovich et al.l (2001). 
The limiting surface brightness of rest-frame optical im- 
ages of LBGs in the NICMOS HDFN indicates that these 
features are likely often too faint to be detected. Also, 
the rest-frame wavelength of the reddest NICMOS im- 
ages of LBGs (i/i6o-band) corresponds to only about 
half that of our z-band observations of the UVLGs and 
so are less sensitive to detection of an older population. 
Resolution effects could play a role as well. In paper II 
we will carry out redshift simulations of a large set of 
rest-frame optical images of UVLGs to investigate these 
optical structures in detail. 

Because high surface brightness regions in LBGs are 
very likely to be blends of super starburst regions, 
we stress that one must be cautious in interpreting 
size/morphology measurements in terms of disk- or 
bulge- like components. The half-light or effective ra- 
dius may be more directly related to the typical radius 
over which the compact, super starburst regions are dis- 
tributed, rather than to the scale size of a (forming) 
bulge or disk (see also lNoguchil [l999: Im meli et al.ll2004t 
lLaw et aT1l2007bf ) , although the two are likely closely re- 
lated if stars predominantly form inside these burst re- 
gions and relaxation or tidal disruption spreads them out 
over time. Although we do not rule out the possibility 
that some fraction of LBGs has bulges at redshifts as high 
as z ~ 4, the high physical resolution data presented here 
suggests that this is possibly far less common than cur- 
rently believed based on morphological parameters (e.g. 
Lotz et al. 2006 estimate -30%). 

Evidence that the small-scale structure of LBG knots is 
indeed of the nature advocated in this paper is found in a 
few rare cases where LBGs are observed at high magnifi- 
cation due to gravit ational lensing. By studying a lensed 
galaxy at z = 4.92, iFranx et all (|1997[ ) were the first to 
show that high surface brightness super starburst regions 
having sizes of a few hundred pc form an essential contri- 
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bution (—75%) to the UV flux in LBGs. Similar features 
can be seen in other lens ed systems that have been dis- 
covered at high redshift (iBunker et alJl200Ct lEllis et al.1 
[200ltlSmail et al.ll2007t [Swinbank et al.ll2007h . We note 
that some fraction of LBGs with very knotty structures 
may be discarded in current large surveys, because their 
light profiles are similar to those of stars. Our sample in- 
dicates that this number could be on the order of 13-25% 
(1-2 out of 8 objects). If we compare the substructure 
of the super starburst regions in Haro 11 and VV 114 at 
z = 0.02 with the same regions simulated at z = 0.15 
(Fig. [6|), we see that majority of the bright nuclei seen in 
the images simulated at z — 0.15 arc a blend of smaller 
super star clusters identified in their unredshifted im- 
ages. However, the most Northern and brightest knot 
in the redshifted image of Haro 11 is still largely un- 
resolved in its unredshifted image, indicating that sin- 
gle, luminous star clusters are capable of dominating the 
rest-frame UV/optical morphology of LBG-like galaxies 
at virtually an y redshift. 

As shown bv lDahlen et al.l (|2007f ) . the fraction of star- 
forming galaxies having a bulge-like morphology at rest- 
frame 2800A decreases from -30% at z w 2 to -10% 
at z 0.5, illustrating the relative importance of con- 
centrated st ar formation at h igh redshift compared to 
low redshift. iZirm et ail (|2007t ) have pointed out the ex- 
istence of a population of compact, massive galaxies at 
z — 2.5 having effective radii and high stellar mass sur- 
face densities. It is possible that these are the result 
of massive and concentrated starbursts in highly dissi- 
pative, gas- rich mergers at high redshift (z — 3 — 6). 
Extrapolation of our results to higher redshifts suggests 
that very dense stellar 'cores' are indeed actively being 
formed in the LBG population. 

5.2.2. Prediction for the faint end slope of the LBG 
luminosity function 

The UV luminosity function (LF) has now been evalu- 
ated over the entire redshift range from z = — 6. Inter- 
estingly, the faint end slope of the LF has been found to 
flatten with dec reasing redshift from a — 1. 74 at z = 6 to 
- 1.2 at z = (lYan fc Windhorstll2004a|bl: iWvder et all 
120051 iBouwens et al.1 l2007t iRvan et al.ri2007h . Vari- 
ous simulations and models have been tried to re- 
produce the LF and explain its flattening at faint 
magnitudes (e.g. iNight et alj l2006t iFinlator et ail 120071 : 
iKhochfar et"aTTl2007h . As we will show, the exis- 
tence of compact starbursts in LBGs may offer a nat- 
ural explanation for the steep slopes that are observed 
at high redshifts. The LF of star-forming regions is 
well-known to have a slope of a « 2 over a wide 
range in physical scales from HII reg i ons to SSCs 
and b eyond (e.g. iKennicutt et all 119891: TMeurer et al.l 
19951: lEhnegreen fc Efremovi Il997p Zepf et al I 119991 : 
Alonso-Herrero et alj|2002t iBradlev et al.ll2006T and ref- 
erences therein). 

We consider the following, very simple, toy model. Let 
N(L, z) be the total number density of galaxies with UV 
luminosity between L and L + dL and redshift between z 
and z + dz. We will also assume that the stellar IMF, as 
well as the star formation efficiency in galaxies, does not 
change with redshift. The UV luminosity of an individual 
galaxy at a given time t is the sum over the luminosity of 



all its previous star- formation events, Li, at that time: 
£lJg(*)=EW W 

i 

The luminosity arising from an 'event' can take any arbi- 
trary form (e.g. burst, constant, declining, etc.), and it is 
not important what caused the event. For LBGs, it is in- 
structive to make a distinction between the contribution 
to the total luminosity from a series of N (semi-)discrete 
starbursts, -k Durst , and that due to a diffuse star- forming 
component that is slowly evolving with time, -^diffuse, so 
that 

JV 

£lBg(*) = -^diffuse (t) + E L l hurst (t). (5) 

i=l 

We have learned that the UV luminosity from an LBG 
is typically dominated by a series of iVknot 'knots' (—1-2 
kpc in size), and that each knot can be further resolved 
into a series of -/Vssb super starburst regions (—100-300 
pc in size), and thus 

Nknot NSSB 

LY% G {t) = L diSusc (t) + E E L ssb(^ (6) 

* 3 

where Lg^it) is the luminosity of a SSB region j in knot 
i at time t. For (proto-) galaxies that are just undergoing 
their first episode of massive star formation the above ex- 
pression will simplify enormously, as we can neglect the 
diffuse component, Ldiffuse (t) , and both N^ not and Vssb 
will be near unity. The most simple case is represented 
by iVknot = Vssb = 1, as observed in some of our lo- 
cal LBG analogs, and (at least in this extreme case) the 
total UV LF will be purely determined by a random sam- 
pling of the LF of the SSBs over the gala xy population 
in a g iven redshift interval. As shown by iMeurer et al.l 
(1995), the 'anatomy of starbursts' typically consists of 
irregularly shaped cloud(s) (SSBs) with embedded com- 
pact sources (SSCs), the brightest of which are massive 
forming (globular) clusters. Star formation thus occurs 
in hierarchically clustered systems ranging from sub-pc 
to multi-kpc scales, and is believed to be generated at 
all scales by self-gravity and turbulence within an un- 
derlying h ierarchical gas mass distribution n (M)dM oc 
M- a dM (jElmegreen fc Efremovi 1 19961 119971 and refer- 
ences therein). Because the UV luminosity of starbursts 
is dominated by its OB stars at all scales, on average 
M oc L with a — 1.7 — 2 for the references given above. 
Therefore, in the regime where the LF is dominated by 
young galaxies generating their first significant amount 
of stellar mass, i.e. at relatively faint magnitudes and 
high redshift, our simple model predicts that the galaxy 
LF observed will essentially be a random subsampling of 
the LF of starburst regions thereby maintaining its in- 
trinsic slope. It is interesting to note that, for example, 
the faintest -B435, Vgo6! or «775-dropout galaxies observed 
at z — 4 — 6 have stellar masses very similar to those of 
our starburst regions (— 10 8 ~ 9 Mq). 

Generally, it will be impossible to model the full UV 
LF quantitatively in this way (e.g. Eq. [3]) without re- 
sorting to semi-analytic models or simulations because it 
involves the convolution of large numbers of galaxies with 
a range of star formation histories including effects such 
as merging, feedback and dust. However, qualitatively 
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one expects that, when all processes are considered, the 
slope will flatten from its initial value. Thus, we expect 
that the faint end slope of the LF at least steepens with 
rcdshift and will deviate less and less from the 'pure' star 
formation LF with a slope of a ks 2, exactly as observed. 

6. CONCLUSIONS 

The study of galaxies in the early Universe can benefit 
enormously from the study of objects that are relatively 
nearby, provided that suitable samples of local analogs 
can be found. H05 and H07 selected "supercompact" , 
UV-luminous galaxies in a large GALEX-SDSS cross- 
matched sample (see §1). Because these objects match 
LBGs in terms of size, SFR, surface brightness, mass, 
metallicity, kinematics, dust, and UV-optical color, this 
sample of local LBG analogs is well suited for investigat- 
ing the connection between star formation and morphol- 
ogy at a level of detail and precision that is impossible 
for high redshift galaxies. Emission line diagnostics in- 
dicate that they are dominated by starbursts, and that 
the line ratios show small offsets with respect to the local 
star-forming population analogous to those observed for 
high redshift LBGs (§3.2). 

In this paper we present HST imaging data in r/330, 
Z850 and Ha of a subsample of 8 local LBG analogs in or- 
der to investigate their morphologies. The effective radii 
estimated from the U330 images range from unresolved 
(<100 pc) to ~ 2 kpc, confirming that the objects are 
highly compact. The objects are slightly more extended 
in 2850 (~l-2 kpc), and they generally have color gra- 
dients with relatively blue inner colors and redder outer 
colors. These gradients are due to the fact that the star- 
burst regions are very compact and are embedded in an 
older, more extended population (§§3.3-3.4). The inter- 
nal extinctions amount to E(B — V) « 0.01 — 0.14 mag 
based on the bolometric dust to FUV luminosity ratios 
observed with Spitzer and GALEX (§3.5). We have cal- 
culated the morphological parameters G, M20 and C and 
find that local LBG analogs have irregular morphologies 
that, on average, lie in between those expected for major 
mergers, disks and bulges. We have simulated our data at 
the depth and resolution of LBGs at z = [1.5, 3.0, 4.0] in 
the COSMOS, GOODS, and UDF surveys, and we have 
carried out direct comparisons with the morphologies of 
z ~ 1.5 star-forming galaxies and z ~ 4 LBGs extracted 
from these surveys. The morphological parameters mea- 
sured for the artificially redshifted sample are very sim- 
ilar to those obtained for the distant LBGs. This estab- 
lishes that local LBG analogs and high-redshift LBGs are 
likely to be in a similar phase of their evolution (§3.6). 
The UVLGs are more typical of LBGs in comparison with 
two very nearby LBG analogs, Haro 11 and VV 114 at 
z = 0.02. The latter two are more similar to the "clump 
cluster galaxies" due to their multiple n uclei at relatively 
large separation (jElmegreen et al.ll2007f ). Although more 
local LBG analogs of this type may be found as well, we 
note that our current sample is biased against finding 
such objects (§3.6.4). 

Although a large fraction of LBGs at z ~ 3 — 4 is 
consistent with having exponential or r 1 / 4 light pro- 
files, their visual morphologies include multiple, dou- 
ble and perturbed systems. Morphological parameters 
indicate that the majority of LBGs are likely to be 
in a range of stages associated with minor and major 



galaxy merging (e.g. iGiayalisco et al.lll996afc lLotz et al.l 
120061 [Ravindr anath et all 120061: lElmeereen et al J 120011 



Burgarella et al.l (|2006h nnd that the ma- 



in contrast, 

jority of luminous, UV-selected star-forming galaxies at 
z ~ 1 has a disk or spiral morphology. A t a similar red- 
shift, but at much lower UV luminosities. Ide Mello et aTl 
(|200rl find a mix of morphological types that include 
early-type galaxies as well as starbursts. We note that 
this is a direct consequence of the selection of galaxies 
purely based on their UV emission. As discussed here 
and by H05, at relatively low redshifts (z ~ — 1) the 
construction of a sample of galaxies that most closely 
resembles high redshift LBGs in terms of their physical 
properties requires additional selection criteria (e.g. UV 
surface brightness) rather than selection on UV luminos- 
ity or the presence of a 'Lyman break' alone. 

We have carried out a detailed investigation of the 
small-scale structure of the local LBG analogs in terms 
of their star forming regions (§4). The starburst regions 
observed in the local LBG analogs are highly compact 
(radii of 100-300 pc), and characterized by bright, unre- 
solved knots of emission within a larger region of diffuse 
star formation that extend up to a few kpc in size. The 
total UV emission from these compact super starburst 
regions (SSBs) is substantial, indicating that 30-80% of 
the total SFR is generated in these knots. Comparison 
with STARBURST99 evolutionary tracks indicates that 
some of the regions are due to very young (~6 Myr) and 
massive (~ 10 s Mq) bursts, while other regions may have 
been forming stars in a more continuous manner for sev- 
eral tens to several hundreds of Myr as indicated by their 
relatively red colors and high Ha equivalent widths. The 
super starburst regions are likely to be a blend of dif- 
fuse star forming regions, stellar associations and super 
or globular star clusters, as demonstrated by comparing 
the images of Haro 11 and VV 114 at z = 0.02 to the 
same images redshifted to z = 0.15 — 4.0. The SSBs are 
generally embedded in a diffuse, older component having 
ages ranging from a few Gyr to a Hubble time, depend- 
ing on the star formation history. In this respect, the 
precursors of the UVLGs may have been more evolved 
than those of LBGs at high redshift (although the pres- 
ence of a "maximally old" stellar component in LBGs is 
consistent with the observations). 

Although some of the SDSS images showed limited ev- 
idence for close companions or faint extended emission, 
the optical HST images clearly reveal disturbed mor- 
phologies at scales well below the SDSS seeing or sensi- 
tivity. Each object shows evidence of merging or interac- 
tions, such as multiplicity of position angles, tidal debris 
or tails, or close companions. Most morphological in- 
formation about high redshift galaxies has been derived 
from rest-frame UV images. In our sample, the E/330- 
band does not establish a definite case for a merger in any 
single one of our galaxies. The disturbed optical (z-band) 
morphologies, together with the luminous, compact su- 
per starburt regions, suggests that mergers of relatively 
gas-rich objects trigger vigorous episodes of star forma- 
tion, with general properties reminiscent of simulations 
of collisional starbursts (§5.1). Our results on local LBGs 
and their similarity to their high-z LBGs constitute the 
most direct evidence to date that the onset of star forma- 
tion and morphological structures of high redshift LBGs 
are largely driven by highly dissipational merging. 
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We discuss several implications of our results for the 
interpretation of lower resolution data on high redshift 
LBGs. Resolution effects will cause the compact star- 
burst regions to blend into typically one or two bright 
LBG-like 'knots'. Although these knots will appear rel- 
atively featureless with total half-light radii of 1-2 kpc, 
the substructure of our LBG analogs indicates that they 
should generally not be interpreted as being evidence of 
bulges or disks. Evidence that the substructure of high 
redshift LBGs is indeed dominated by SSBs is found in 
a few cases where LBGs are highly magnified by grav- 
itational lensing (§5.2.1). However, the strong relation 
we find between merging and the triggering of very com- 
pact, massive super starburst regions indicates that these 
events may be closely linked to the formation of stellar 
bulges. Furthermore, we suggest that the prominence of 
luminous, unresolved (super) starburst regions in form- 
ing galaxies may provide a natural explanation for the 
value of the faint end slope of the UV luminosity function 
at high redshift by an extension of the local star-forming 
region luminosity function which is well-known to have 
a power law slope of a « 2 (§5.2.2). 

Future work. - In Paper II we will study a significantly 
larger sample of local LBG analogs to be observed with 
HST in Cycle 16 in the FUV/optical using the ACS/SBC 
and WFPC2. Follow-up spectroscopy is being used to 
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TABLE 1 

Sample and log of HST/ACS observations. 
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TABLE 2 

Star formation rates and stellar masses. 
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7.3 


<3.7 


2.14 


SDSS 


J092600.41+442736.1 


6.8 ±0.4 


3.7 


7.8±1.4 


0.16 


SDSS 


J102613.97+484458.9 


3.4 ±0.1 


6.6 


9.2±1.3 


0.81 


SDSS 


J135355.90+664800.5 


7.1 ±0.9 


14.6 


22.5±2.7 


0.98 


SDSS 


J214500.25+011157.6 


3.4 ±0.8 


7.5 


14.5±2.9 


1.32 



TABLE 3 

Photometry and UV spectral slope. 



Galaxy 



mpuv 
(mag) 
19.23±1',^ 
-,+0'.2i 
-0.31 
(-0.06 
J -0.06 
f0.10 
-0.11 
,+0.11 
3 -0.13 

19.271°;°! 
18.99lS-.8J 
19-94l°;i 



m N uv 
(mag) 

18.85lH« 
19.231°;^ 



m,33o m u rosso Pfuv-nuv 
(mag) (mag) (mag) 



SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 



J005527. 
J032845. 
J040208. 
J080844. 
J092600. 
J102613. 
J135355. 
J214500. 



46-002148.7 
99±011150.8 
86-050642.0 
26±394852.4 
41±442736.1 
97+484458.9 
90±664800.5 
25±011157.6 



19.42] 
18.90j 
18.22^ 
18.83^ 



hO.03 
J -0.03 
hO.05 
-0.05 
3+0.07 
-0.08 

19.00l°;° 2 2 

7+0.02 



18.70] 
17.59j 
18.89"* 



18.57_ 



0.02 



19.22±g;i§ 



18.78 
18.97 
18.81 
17.62 

_a 

18.86 
18.42 
19.61 6 



18.65 ±0.02 
18.68 ±0.04 
18.67 ±0.03 
17.48 ±0.01 
19.04 ±0.03 
18.85 ±0.02 
18.35 ±0.02 
19.24 ±0.04 



18.25 
18.15 
18.08 
17.05 
19.28 
18.23 
17.67 
18.62 c 



-1.13 ±0.22 
-1.56 ±0.39 
-1.54 ±0.09 
-0.54 ±0.15 
-2.12 ±0.18 
-1.38 ±0.05 
-1.03 ±0.06 
-0.35 ±0.41 
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TABLE 4 
Size measurements and colors 



Galaxy 


r 5 o,M 


^50,330 


r 90,330 


r 50,850 


r 90,850 


r50,Ha 




(f330-Z85o)in ne 


(E/330-2850)outer 




(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(mag) 


(mag) 


SDSS ,1005527.46-002148.7 


0.77 


0.37 


1.26 


0.78 


2.13 


0.50 


1.69 


0.03 


0.89 


SDSS J032845. 99+011150. 8 


1.66 


0.87 


1.96 


1.79 


4.93 


0.93 


2.49 


0.24 


1.32 


SDSS J040208. 86-050642.0 


1.16 


0.84 


1.94 


1.44 


4.26 


1.23 


2.81 


0.48 


0.84 


SDSS J080844. 26+394852. 4 


0.47 


0.16 


1.41 


0.57 


2.77 


0.47 


2.94 


-0.58 


1.52 


SDSS J092600. 41+442736.1 


0.97 






1.10 


3.33 


0.87 


2.67 






SDSS J102613. 97+484458. 9 


2.07 


1.89 


3.86 


1.99 


4.48 


2.13 


4.22 


0.53 


0.67 


SDSS J135355. 90+664800. 5 


1.81 


1.46 


3.19 


3.56 


9.06 


2.74 


7.23 


0.48 


0.81 


SDSS ,1214500.25+011157.6 


0.74 


1.10 6 


3.55 6 


1.30 c 


4.64 c 


d 


d 


0.93 


1.00 



TABLE 5 

Infrared photometry and extinction. 





Galaxy 


/3.6 


h 




/24 


/70 


/16O 


TI R 


n FUV 


Ad 
^FUV 


e(b - vy int 








(pJy) 




(mjy) 


(mjy) 


(mjy) 


(10 10 L ) 


(10 10 L ) 


(mag) 


(mag) 


SDSS 


J005527. 46-002148. 7 


725.8+33.2 


4280.5+80.8 


33 


2+0.4 


163.4+6.4 


115.4 


16 


3.2 


1.48 


0.14 


SDSS 


J032845. 99+011150. 8 


218.3+39.0 


972.0+71.9 


3 


4+0.4 


<23.4 / 


25.9 


<1.9 


1.8 


0.56 


0.05 


SDSS 


J040208. 86-050642.0 


131.6+34.8 


415.3+66.6 


1 


1+0.1 


<19.S f 


25.5 


<1.6 


2.8 


0.31 


0.03 


SDSS 


J080844. 26+394852. 4 


1048.6+39.8 


5786.5+66.2 


39 


9+0.6 


156.0+6.8 


117.0 


4.1 


2.0 


0.91 


0.08 


SDSS 


J092600. 41+442736.1 


90.7+4.2 


229.4+24.0 


3 


4+0.4 


20.7+3.7 


10.2 


2.1 


5.5 


0.16 


0.01 


SDSS 


J102613. 97+484458. 9 


176.6+36.2 


614.5+56.8 


4 


2+0.3 


42.5+4.2 


31.6 


3.7 


2.7 


0.70 


0.06 


SDSS 


J135355. 90+664800. 5 


369.2+32.4 


1660.3+53.1 


44 


4+0.4 


70.1+4.3 


27.2 


8.2 


6.0 


0.71 


0.07 


SDSS 


J214500.25+011157.6 


200.4+31.5 


1617.5+64.2 


5 


0+0.4 


23.0+5.5 


27.7 


4.2 


2.6 


0.79 


0.07 



TABLE 6 
Morphological parameters 





Galaxy 


^P,330 
(") 


G330 


-^20,330 


C330 


rp,850 
(") 


^850 


^20,850 


C850 


SDSS 


J005527.46-002148.7 


0.55 


0.63 


-2.41 


4.73 


0.83 


0.58 


-1.98 


3.53 


SDSS 


J032845.99+011150.8 


1.66 


0.66 


-1.48 


3.97 


2.62 


0.58 


-1.55 


3.34 


SDSS 


J040208.86-050642.0 


1.02 


0.56 


-1.60 


3.29 


1.92 


0.60 


-1.98 


3.92 


SDSS 


J080844.26+394852.4 


0.13 


>0.57 


<-2.55 


>3.01 


2.97 


0.66 


-3.23 


7.11 


SDSS 


J092600.41+442736.1 


_a 


_a 


_a 


_a 


1.29 


0.73 


-1.60 


3.27 


SDSS 


J102613.97+484458.9 


2.21 


0.60 


-1.31 


3.28 


1.92 


0.58 


-1.48 


2.93 


SDSS 


J135355.90+664800.5 


1.26 


0.60 


-1.27 


2.88 


1.29 


0.53 


-1.57 


2.88 


SDSS 


J214500.25+011157.6 


0.75 ft 


0.56 6 


-1.47 6 


2.95 ft 


1.40 c 


0.61 c 


-1.80 c 


3.56 c 
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Fig. 1. — Panels show Galaxy Evolution Explorer (GALEX) in the far-UV (FUV) and near-UV (NUV) of the 8 objects in our sample. 
Bars in the lower right of each panel measure 5", and are comparable to the full- width-at- half-maximum of the GALEX point spread 
function. North is up, East is to the left. 






Fig. 2. — Panels indicate the HST morphologies in the UV (left), optical (middle), and Ho (right). The images were smoothed using a 
Gaussian filter of 1.5 pixels (FWHM) to reduce some of the noise at the sub-PSF level, and shown using a logarithmic stretch. The physical 
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Fig. 3. — Optical emission line diagnostic diagrams log([OIII]/H/3) vs. log([NII]/Ha) (left panel), log([OI]/Ho) (middle panel), and 
log([SII]/Ha) (right panel). Points indicate all objects in the GALEX Medium Imaging Survey (MIS) Galex Release 2 (GR2) sample 
cross-matched with the SDSS DR4 sample. The objects that are the subject of the current paper are indicated by the large filled circles. 
Small filled circles indicates the other objects in our main sample of compact UVLGs that have HST, Spitzer and radio data either taken 
or sch eduled. The solid line indicates the boundary between regions typically populated by starbursts and Seyferts from Kauffman n et alj 
(2003). Dashed lines indicates the boundary between objects of starburst spectral type and pure AGN from lKewlev et al.l 112009) . Objects 
between the solid and dashed lines in the left panel have typically composite spectra consisting of a metal rich stellar population and 
an AGN. Note that, based on these diagnostics, object J080844. 26+394852. 3 is peculiar. The standard diagnostics diagram (left panel) 
indicates it could be of composite type, while it lies well on the sta rburst side in the o ther diagrams (middle and right panels). This 
behavior is typical for galaxies with very high ionization parameters (Kcwley et al. 2001). The remainder of our HST objects lies in the 
starburs t dominated region, albeit offset towar ds relatively high values of log([OIII]/H/3. A similar offset has been found for LBGs at high 
redshift l|Shaplev et al. 2005; Erb et al. 200H|). 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

Radius (arcsec) 

Fig. 4. — UV-optical radial color profiles of the UVLGs shown in Fig. [2] Lines indicate the (cumulative) color determined in circular 
apertures measured from the centroids of the object detected in the optical image out to the 90% flux radius (dashed colored lines), and with 
respect to the object centroids defined in the UV image (dotted colored lines). The UVLGs are generally very blue within approximately 
the half light radius, and become redder at larger radii. The objects often have steeper inner color gradients when then UV image is used 
for object detection (dotted lines), compared to when the optical image is used (dashed lines) due to the different positions of the object 
centroids in the UV and optical. For 214500, the UV-optical color corresponds to FUVi50-V@o6, whereas for the other objects the colors 
measured are C/330-2850 • Errors are 3<r. 
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Fig. 5. — Distribution of morphological parameters M20 vs. G (left panel) and G vs. C (right panel) for the sample of z ~ 0.15 UVLGs. 
Morphologies measured from the UV images are indicated by stars, morphologies measured from the optical images by filled circles. Dotted 
lines connect th e UV and optical morphologies measured for each source. For comparison, we have indicated the following (all adapted 
from ILotz etaD (2006)): regions populated by likely mergers having multiple nuclei (left hatched area) and by bulge-dominated objects 
(right hatched area), regions occupied by simulated de Vaucouleurs and exponential profiles (top and bottom boxes, resp.), the median 
values found for star-forming galaxies at z ~ 1.5 and LBGs at z ~ 4 in GOODS (encircled crosses), and the typical errors for the high 
redshift samples (large crosses). 
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Fig. 6. — Top panels: NUV archival images of the local starburst galaxies Haro 11 and VV 114, both at z = 0.02. The image of Haro 
11 was taken with the ACS/HRC (Program 10575, PI: Goran Ostlin). The image of VV 114 was taken with STIS (Program 8201, PI: 
Gerhardt Meurer). See Sect. l2.4.1l for details. Bottom panels: Simulated F330 W ima ges of Haro 11 (left) and VV 114 (right) at z = 0.15 
for comparison with our UVLG sample. The images measure 3"x3". See Sect. I3.6.2l for details. 



New evidence for a merger-driven formation of LBGs 



25 



005527 



* 

F330W 



0328J5 



F330W 



04BS0S 
F330W 



• 

F330W 



102613 



FJ30W 



13535.5 



~33t\Y 



214500 
F150LP 




Kara 11 




= 3.0 



is « 3.0 




Fig. 7. — Redshift simulations of the UVLG sample at the depth of GOODS (Sect. 13. 6. 21 . Panels (3" X 3") show the original, unredshifted 
HRC/F330W or SBC/F150LP images (left), and GOODS B 435 or V 6 q 6+»775 simulations at z = 1.5, z = 3.0 and z = 4.0. The last two 
rows show the redshift simulations of Haro 11 and VV 114 (Sect. 13. 6.21 1 for comparison. 
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Fig. 8. — Rcdshift simulations at the depth of the UDF. See Sect. 13.6.21 and the caption of Fig. [7]for details. 
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Rcdshift simulations at the depth of COSMOS. See Sect. I3.6.2l and the caption of Fig. [7]for details. 



28 



Overzier et al. 




-3.0 1 1 1 1 . ' * 




-0.5 h i i i 

12 3 4 

Redshift (z) 




Fig. 10. — Morphological parameters of the UVLGs measured from the simulated images at different redshift. Objects were redshifted 
to z = 1.5, 2 = 3 and 2 = 4, and simulated at the depths of COSMOS (red dashed lines), GOODS (green dotted lines) and the UDF (blue 
solid lines). 
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Fig. 11 . — Co mparison of the morphologies of Lyman break galaxies at ~ 4 and starburst galaxies at z ~ 1.5 from the samples of 
ILotz et al.1 (120061 ) (plusses) with the morphologies of the UVLGs (stars). The UVLGs were artificially redshifted to the mean redshift of 
the comparison samples, and simulated at a similar depth. Top panels show the results for LBGs and UVLGs at the UDF depth at z ~ 4, 
middle panels show LBGs and UVLGs at the GOODS depth at z ~ 4, and bottom panels show star forming galaxies and UVLGs at the 
GOODS depth at z ~ 1.5. Haro 11 ('H') and VV 11 4 ('V') wer e simul ated at the same redshift and depth as the UVLGs. For comparison, 
we have indicated the following (all adapted from lLotz et all ^2006ft ): regions populated by likely mergers having multiple nuclei (left 
hatched area) and by bulge-dominated objects (right hatched area), and the typical errors for the high redshift samples (large crosses). 
In all panels, the UVLGs span a very similar range in parameter space compared to the comparison samples, indicating that UVLGs and 
LBGs have very similar quantitative morphologies. 
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Fig. 12. — Panels show the C330 images indicating the apertures placed around the unresolved star burst regions identified by eye in each 
of the sources (small circles). Apertures measure 0'.'3 in diameter, and are numbered for identification of each region in Fig. 1141 Larger 
apertures were chosen to measure the average color of the regions in which the unresolved starbursts are situated (large circles). For object 
214500 the image shown is FUVibo- F° r completeness, we show the 2:850 image for 092600 as no C/330 image is available. 
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Fig. 13. — Apparent (left panel) and absolu te (right panel) magnitude distributions in C/330 (solid lines) and 2850 (dotted lines) of the 
unresolved starburst regions identified in Fig. 1121 The magnitudes were extinction-corrected using a global reddening value derived from 
the bolometric dust to FUV luminosity ratio (see Table [3). 
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Fig. 14. — Color magnitude diagram of the individual starburst regions identified in Fig. 1121 Filled, numbered symbols correspond to the 
unresolved regions identified in 005527 (circle), 032845 (star), 040208 (upward triangle), 080844 (downward triangle), 102613 (ellipse), and 
135355 (square). The large, open symbols correspond to the large apertures in Fig. I12l with the point sources removed. Tracks indicate the 
color and magnitude evolution of an instantaneous burst model (left panel) and a continuous star formation model (right panel) simulated 
using STARBURST99. Blue tracks are for LMC metallicity (Z = 0.008), red tracks for solar metallicity (Z = 0.02). For instantaneous 
bursts we plot tracks having a burst mass of 10 8 Mq and 10 7 Mq. For continuous models we plot tracks having SFRs of 5, 1, and 0.1 
Mq yr -1 . Ages in Myr have been indicated along the tracks (marked at 6,7,8,9,10,12,14,16,18,20,30,40,50,60,70,80,90,100 Myr for the 
burst models, and at 10,20,30,40,50,60,70,80,90,100,200,500,1000 Myr for the continuous models). Dotted lines indicate the colors if we 
include the contribution from the nebular continuum. The model tracks are reddening-free, but the observed colors and magnitudes were 
dereddened using a global reddening value derived from the bolometric dust to FUV luminosity ratio (see Table [3j . The starburst regions 
of 214500 were omitted because of the different filterset used during those observations. Object 092600 was omitted because no r/330 is 
available. 
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Fig. 15. — Rest-frame equivalent width Ha vs. C/330— £850 color for the super starburst regions. Tracks indicate the age evolution of the 
Ha EW for instantaneous burst (dashed lines, with ages marked at 5,6,7,8,9,10 Myr) and continuous star formation (solid lines, with ages 
marged at 5,6,7,8,9,10,20,30,40,50,60,70,80,90,100,200,500,1000 Myr) models from STARBURST99. See text and the caption of Fig. [Til for 
further details. 
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APPENDIX 
SCARPA ET AL. (2007) REVISITED 

In a recent paper, IScarpa et al.l ( 2007) h ave tried to ar gue that the sample of compact UVLGs selected by H05 and 
H07 are not true analogs of LBGs. IScarpa et al.l ([2007) claim that the sizes estimated from the SDSS u-band images 
grossly underestimate the actual size in the UV, and that the UVLGs therefore are much larger and lower in surface 
brightness than LBGs. Here, we will briefly show that this claim is based on the combination of a faulty analysis and 
incorrect assumptions, that are easily refuted by the HST data presented in this paper. 

• IScarpa et al.l (|2007l ) claim that it is only possible to detect central bulges in the SDSS M-band images, and that 
fainter, extended disks would be missed. They reach this conclusion by fitting the SDSS r-band images with a 
two component model (point source plus disk) and then using the half-light radius of only the disk component to 
recompute the size of the UVLGs. In doing this, the authors made a serious mistake, because they ignore the fact 
that on average, half the light in their model fits comes from the central "point source" (see their Table 1). 

• The claim in the IScarpa et al.l (|2007f) paper that most of the FUV flux in compact UVLGs should come from a 
large disk component relies on the observation that there exists a population of galaxies with relativ ely red bulges 
not visible in the FUV and very blue outer rings or disks that are only apparent in the UV (e.g. iThilker et alJ 

2005) . It should be clear both from the results presented in H05 and H07, and from the analysis of HST data 
presented here, that UVLGs are an entirely different class of objects. In fact, the UVLGs are clearly more 
extended in the optical than in the UV due to an underlying (older) component. The UVLGs have very blue inne r 
colors, and become redder outwards (see Fig. [4]and TablelJ)- This is opposite to the assumption of lScarpa et al.l (|2007h . 

• The HST C/330 images confirm that the UVLGs are indeed very compact objects. In fact, in most cases, the 
half light radius is even smaller than the seeing deconvolved SDSS u-band size (Table |4|). In the case of both 
J005527.46-002148.7 and J080844.26+394852.3 the UV flux is dominated by an unresolved component even at the 
HST resolution. Furthermore, when we extrapolate the slope of the UV continuum measured from the FUV-NUV 
color and predict the 6^330 flux, we find good agreement with the measured 6^330 flux (Sect. 13. 3j) . indicating that the 
£^330 (and SDSS w-band) image is a good tracer of the FUV morphology. In the case of J214500. 25+011 157.3 we have 
an A CS/SBC image taken through F150LP, a filter that is almost identical to the GALEX FUV filter (see lTeplitz et al.l 

2006) . The corresponding FUV surface brightness calculated solely from the ACS data, log w lFUV,ACS = 9.53 Lq 
kpc -2 , is slightly lower than the value of H07 based on GALEX and SDSS (\og 10 Ipuv = 9.82 L© kpc~ 2 ), but still 
well above our LBG selection threshold of \og 10 lFuv = 9 Lq kpc -2 . 

• As we have discussed above, the galaxy sizes defined by IScarpa et al.l (|2007f ) pertain only t o the disk component and 
thus cannot be compared with the half-light radii quoted by H05 and H07. IScarpa et al.l (|2007l ) also did not apply 
their profile fitting methodology to the full sample in a uniform way. They only applied their analysis to the objects 
that were marked as "supercompact" by H07 and not to the other UVLGs. As a result, their concluding Figure 5 is 
misleading. It appears to show that the LBG analogs have a similar FUV surface brightness as the main sample. Had 
they applied their technique to all the objects in the sample, the compact UVLGs would still form an upper envelope 
to the Ifuv distribution, with the overall Ifuv distribution shifted downwards toward lower surface brightness due 
to the fact that the "point source" flux in their fits was ignored. 

• Finally, we would like to remind the reader that the real test of whether or not a given sample of nearby galaxies 
are good analogs of LBG s is not the UV select ion per se (note that the LBG luminosity function extends over many 
orders of magnitude (e.g. iBouwens et"aT1l2007l )). but whether their main physical properties are indeed similar. We 
firmly believe that our sample of LBG analogs have passed this test. 



